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I. INTRODUCTION 



A . BACKGROUND 

One of today's major environmental concerns is the depletion of the 
earth's protective ozone layer. In 1987, an international conference was 
held in Montreal, Canada to address the problems caused by 
Chlorof luorocarbons (CFCs) to the earth's ozone layer. CFCs are manmade 
chemicals of chlorine, fluorine, and carbon and are unique in that they 
have a combination of desirable properties: low in toxicity, non- 
flammable, non-corrosive, non-explosive, extremely stable and compatible 
with many other materials. This extreme stability is what causes problems 
to the ozone layer due to the fact that CFCs only break down in the upper 
atmosphere when subjected to intense ultraviolet radiation. This break 
down produces chlorine which has been linked to the depletion of the 
earth's ozone layer. 

In September 1987, 24 nations representing the United Nations 
Environment Program (UNEP) met and signed the Montreal Protocol. They 
discussed the substances that deplete the ozone layer [Ref.l] and called 
for a near-term freeze on the production and consumption of these 
substances. The agreement required production of these chemicals to be 
cut back to 1986 levels followed by a two-phased reduction culminating in 
cutbacks of 50% by mid-1998; this came into effect on July 1, 1989. In 
1990, a progress meeting was held in London where UNEP delegates agreed to 
completely phase out all CFCs by the year 2000 [Ref. 2]. In the spring of 
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1992, President Bush pushed up the complete phase out od CFC's by the year 
1995. 

The U.S. Navy uses a number of different CFCs (designated by "R" for 
refigerants) for various refrigeration and air conditioning (AC) needs. 
Presently, the U.S. Navy has approximately 1850 shipboard AC plants using 
both R-12 (in reciprocating compressor) and R-114 (in centrifugal 
compressor) plants. To comply with the Montreal Protocol and U.S. 
legislation, the Mechanical Systems Branch (Code 2772) at the Naval 
Surface Warfare Center (NSWC) is pursuing research in the elimination of 
shipboard use of CFCs. As mentioned by Chilman [Ref. 3] this research is 
to be completed in three phases: 

1. To identify in the short term suitable alternative ozone-safe 
chemicals to replace R— 114 and R— 12. To accomplish this task, the heat 
transfer characteristics must be similar to the existing refrigerants in 
place and hence the need for a database exists for current refrigerants 
(R-114 and R-12) so that they can be compared to the new proposed 
refrigerants (HFC-124 and HFC-134A respectively) . 

2. In the longer term, to research, develop, and test substitute 
chemicals and alternative technologies to replace existing CFC uses. 

3. To implement new cooling system technologies into the fleet which 
do not depend on CFCs or their replacement. 

This thesis is a continuation of the previous work at NPS and 
supplements NSWC's research on alternatives to CFCs by establishing 
baseline nucleate pool boiling data of pure R-114 and R-114/oil mixtures 
from a small bundle of enhanced tubes representing a section of a flooded 
evaporator. Emphasis is placed on the natural convection and boiling 
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regions, hystersis phenomena, and analysis of various oil concentrations 
on the overall heat transfer performance. 

B. OBJECTIVES 

The objectives of this thesis are as follows: 

1. Understand in greater detail both the convection and nucleate pool 
boiling phenomena and hysteresis effects within a small Turbo— B tube 
bundle . 

2. Obtain data using a Turbo-B tube bundle for increasing and 
decreasing heat flux for R-114/oil mixture with oil concentrations of 0, 
1, 2, 3, 6, and 10 percent. 

3. Compare data with earlier studies at the Naval Postgraduate School 
using R-114/oil mixtures for other enhanced tube surfaces. 
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II. LITERATURE SURVEY 



A. SMOOTH TUBE BUNDLES 

In recent years significant progress has been made in understanding 
nucleate boiling heat transfer phenomena on the shell side of flooded 
evaporators. Extensive work on smooth tube bundles has been reported by 
Cornwell (Leong and Cornwell [Ref. 4], Cornwell et al. [Ref. 5], Cornwell 
and Scoones [Ref. 6], Cornwell [Ref. 7]). Cornwell and Schuller [Ref. 8] 
conducted a photographic study of boiling R-113 in a smooth tube bundle at 
one atmosphere. One of their conclusions was that bubbles leaving the 
lower tubes in the bundle impacted and caused a sliding motion around the 
upper tubes. Cornwell and Schuller observed the two-phase flow patterns 
and deduced that sliding bubbles from lower tubes on upper tubes could 
account for significant heat transfer in the top part of the bundle. 
Cornwell [Ref. 7] later found that in the nucleate boiling region, sliding 
bubbles and liquid forced convection could account for all the heat 
transfer in the top of the bundle. 

The influence of tube position within a bundle of smooth tubes using 
R-ll, R-12 , R-22 and R-113 has been studied extensively by Wallner [Ref. 
8], Fujita et al . [Ref. 9], Chan and Shourki [Ref. 10], Rebrov et al . 
[Ref. 11], and Marto and Anderson [Ref. 12]. Using both in-line and 
staggered tube arrangements with various tube pitch-to-diameter ratios 
between 1.2 and 2.0, their work verified that the influence of the lower 
tubes in a bundle can significantly increase the heat transfer performance 
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of upper tubes at low heat fluxes due to two-phase convective effects. At 
high heat fluxes (typically > 50 kW/m 2 ) in the fully developed nucleate 
boiling region, the data for all the tubes merged onto a single curve. 
This is representative of a single smooth tube and shows that there is no 
'bundle effect' (ie. no improvement over a single tube under similar 
conditions) in the high heat flux (nucleate boiling) region. 

Chan and Shoukri [Ref. 10] studied the boiling characteristics of a 
smooth in-line tube bundle in R-113. They concluded that at lower heat 
fluxes, the heat transfer process is strongly influenced by two-phase 
convection effects, resulting in higher heat transfer coefficients on the 
upper tubes. At high heat fluxes, however, they found that the dominant 
mode of heat transfer was nucleate boiling from the upper tubes and that 
convective effects from below were insignificant. At these high fluxes, 
the bundle performance was similar to the trends of a single tube in a 
single tube apparatus. Fujita et al . [Ref. 9] also found that the heat 
transfer at low heat fluxes using a smooth tube bundle in R-113 was 
enhanced by convection induced by rising bubbles (ie. a steady increase in 
performance of the upper tubes as additional lower tubes were activated) . 
They attributed this enhancement to the "positive bundle effect". A high 
heat fluxes, this enhancement disappeared. 

Anderson [Ref. 13] found similar effects as above for a smooth tube 
bundle in pure R-114. Furthermore, he reported that the presence of up to 
3% oil (by mass) actually improved the heat transfer performance. This is 
similar to data reported for a single smooth tube by Wanniarachchi et al . 
[Ref. 14]. Furthermore, at an oil concentration of 10%, only a slight 
degradation in the heat transfer (compared to pure R-114) was found. He 
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obtained a maximum heat transfer performance for the bundle at an oil 
concentration of 2% . ^ 

k 

B. ENHANCED TUBE BUNDLES 

Much less work has been done on enhanced tubes (Enhanced means any 
surface that is not smooth). However, the effects at high and low heat 
fluxes mentioned above are similar to those obtained for finned tube 
bundles by Yilmaz and Palen [Ref. 15], Muller [Ref. 16], and Hahne and 
Muller [Ref. 17]. Stephen and Mitrovic [Ref. 18] looked at R-12 and R-114 
boiling from a GEWA-T tube bundle. Apart from the magnitude of the heat 
transfer coefficient varying with fluid, the trends are very similar to 
those mentioned above for smooth and finned tube bundles. 

For porous coated surfaces, Czikk et al . [Ref. 19] found no 'bundle 
effect' over a wide range of heat flux (1-100 kW/m 2 ) and the bundle data 
agreed closely with single tube data. Arai et al . [Ref. 20] found that 
the 'bundle effect' for a Thermoexcel tube bundle was smaller than that 
found for a smooth or finned tube bundle. However as before, any 'bundle 
effect' was eliminated at high heat fluxes where the data for all the 
tubes agreed closely with single Thermoexcel tube results. These effects 
are similar to those found by Czikk et al. [Ref. 19] for the porous 
coated. 

Chilman [Ref. 3] reported experiments with a Turbo-B tube bundle in 
pure R-113, conducting both increasing and decreasing heat flux tests. He 
concluded that in the natural convection region, heated lower tubes do not 
have any influence on the heat transfer from upper tubes. Also, Chilman 
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reported the presence of heated lower tubes within a bundle reduced the 
incipient boiling point. 

Stephan and Mitrovic [Ref. 18] reported the influence of oil on the 
boiling heat-transfer coefficient of R-12 using a T-shaped finned tube 
(Gewa-T) bundle. They reported that the ratio of oil to no oil heat- 
transfer coefficients decreased with the mass fraction of oil for all 
except the highest heat flux (22 kW/m 2 ) where an increase in heat transfer 
was noted for oil concentrations between 1 and 6% . They concluded that 
the influence of oil on heat transfer was mainly due to the thermal 
properties of the specific oil used in the experiments and its interaction 
with the refrigerant. 

Schlager et al . [Ref. 21] summarized the influence of oil on 
refrigerant in pool boiling. They stated that under certain conditions 
(typically low pressure and high heat flux) , the heat transfer coefficient 
increased at low oil concentration. Stephan [Ref. 22] first pointed this 
out and attributed the phenomenon to foaming. Burkhardt and Hahne [Ref. 
23] for a finned tube bundle found that the maximum heat transfer 
coefficient, which was 10% to 15% above the oil-free value, occurred at a 
concentration of about 4%. 

Heimbach [Ref. 24] conducted experiments with R-12/oil mixtures on a 
finned tube bundle. He reported that the presence of up to 2% oil, did 
not affect the heat transfer performance significantly. However at higher 
concentrations (3% to 7%), an increase in the heat transfer was observed. 
He also attributed this to foaming and postulated that changes in the 
properties of the mixture might facilitate the formation of bubbles. 
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Anderson [Ref. 13] and Akcasayar [Ref. 25] conducted experiments with 
finned (19 fpi) and High Flux (porous coated) tube bundles in pure R— 114 
and R-114/oil mixtures in the same apparatus. Anderson reported a maximum 
heat transfer performance at an oil concentration of 3% for the finned 
tube bundle. Akcasayar also reported that the finned tube bundle 
performance increased 1.65 times with 3% oil concentration (compared with 
pure R— 114) at the maximum heat flux level. For 6% and 10% oil 
concentrations, the performance of the bundle, when compared with lower 
oil concentrations, decreased. This was especially significant at a 10% 
oil concentration. Compared with the finned tube bundle, the High Flux 
tube bundle had a 1.5 times better heat transfer performance at a heat 
flux of 30 kW/m 2 in pure R-114. However, these performance ratios 
decreased with increased oil concentrations such that the finned bundle 
outperformed the High Flux bundle at 6% oil concentration. This was 
especially true at the highest heat fluxes where the High Flux bundle 
performance was not much better than a smooth tube bundle of similar size. 
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III. EXPERIMENTAL APPARATUS 



A . TEST APPARATUS OVERVIEW 

The experimental apparatus including the auxiliary equipment and the 
evaporator/condenser is shown in Figure 1. The following is only a 
general description of the whole experimental apparatus. A more detailed 
look at the condenser and evaporator is provided in section C. Further 
information about the apparatus is provided by Murphy [Ref. 26] and 
Anderson [Ref. 13]. 

The experimental apparatus is essentially made up of three closed 
loops. The first loop consists of an 8 ton refrigeration unit located 
outside the laboratory which is used to cool an ethylene glycol/water 
mixture. The second loop is the ethylene glycol/water mixture flowing 
through the condenser. This mixture is contained within a large sump 
within the laboratory. The flow rate through the condenser is delivered 
by two pumps which can be operated independently or together; this coolant 
mixture condenses the refrigerant vapor in the condenser and maintains 
system pressure and temperature. Pump number one provides coolant flow 
through the four test condenser tubes as well as to one of the auxiliary 
condenser coils (bottom coil) . Pump number two provides coolant through 
a manifold which distributes the coolant to the remaining four auxiliary 
condenser coils within the condenser. The third loop is the evaporator 
and condenser itself designed for reflux operation. The vapor generated 
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in the evaporator flows upward and condenses in the condenser; the 
condensate then returns to the evaporator via gravity. 

B. AUXILIARY EQUIPMENT 

1. 28 kW Refrigeration Unit 

This unit was used to cool a 1.8 m 3 reservoir sump of an 
ethylene glycol/water mixture (coolant) to the desired temperature needed 
to condense the refrigerant vapor. For R-114, the temperature control was 
set to maintain the sump at -15 °C . The refrigeration unit had a cooling 
capacity of 28 kW (8 tons). 

2. Ethylene Glycol/Water Mixture 

The coolant used was a 54% (by weight) ethylene glycol/water 
mixture. This refrigerated mixture was used to control the system 
pressure and temperature by circulating coolant through the auxiliary 
condenser coils and/or condenser test tubes at different flow rates. 

3 . Pumps 

Two pumps were available to circulate the coolant from the sump 
through the condenser. Pump number one fed four test condenser tubes and 
one of the auxiliary condenser coils. Pump number two fed the other four 
auxiliary condenser coils. Pump number one was the primary pump used at 
low heat fluxes ; pump number two was started (as necessary) at high heat 
fluxes to maintain the desired saturation pressure in the evaporator. 

4. Flowmeters 

Five calibrated float-type flowmeters, connected to pump number 
one, were used to measure the flow rates passing through the four test 
condenser tubes and one auxiliary condenser coil. One additional 
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flowmeter (connected to pump number two) was used to measure the total 
flow to the remaining four auxiliary condenser coils. Each of the five 
auxiliary condenser coils had a globe valve to regulate (or shut off) 
flow as desired. For the four test condenser tubes, the coolant flow was 
regulated by a flowmeter valve. 

C . EVAPORATOR/CONDENSER 

An overall view of the evaporator and condenser is shown in Figure 2. 
The evaporator was designed to simulate a small portion of a refrigerant 
flooded evaporator. Front and side views of the evaporator are shown in 
Figures 3 and 4. It was fabricated from stainless steel plate and formed 
into a short cylinder, 610 mm in diameter and 241 mm long. Electrically- 
heated tubes were cantilever-mounted from the back wall of the evaporator 
to permit viewing along the axis of the tubes through the lower of two 
viewing windows. A plexiglas plate was attached to the front of the tube 
bundle to ensure tube alignment during experiments. Each viewing window 
had a layer of glass and plexiglas, the glass being used on the 
refrigerant side in order to prevent surface cracking of the plexiglas. 
The plexiglas gave the glass added strength and served as a safety barrier 
in case the glass broke during pressurization. 

The electric power can be applied separately to each set of heaters 
using a STACO 240 V, 23.5 kVA rheostat controller shown in Figure 5. 
Also, the desired number of instrumented tubes, active tubes, simulation 
or auxiliary heaters can be individually activated by using circuit 
breakers. The five simulation heaters, each with a maximum rating of 4 
kW, were mounted below the test bundle in order to simulate 15 additional 
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tube rows in a larger bundle and to provide an inlet vapor quality into 
the bottom of the test bundle as suggested by Webb [Ref. 27] . The liquid 
pool was maintained at 2.2 °C (corresponding to a saturation pressure of 
1 atmosphere) by passing coolant through the condenser. 

Figure 6 is a schematic sectional view of the evaporator that shows 
the four kinds of heated tubes installed in the evaporator. These were 
instrumented, active, auxiliary, and simulation. For this study only the 
instrumented, active, and simulation heaters were used; the auxiliary 
heaters are needed for experiments either at higher pressures or for other 
refrigerants which have a higher normal boiling point (such as R-113) . 
Table 1 gives the power rating for these heaters and the number used in 
the evaporator. 

The tube bundle itself consists of instrumented, active, and dummy 
tubes. The location of each tube is represented by the respective letter 
I, A, and D as shown in Figure 6. The test bundle consists of two types of 
heated tubes: 12 active tubes (marked "A”) which contained 1 kW cartridge 
heaters, and 5 instrumented tubes (marked "I") which, in addition to the 
1 kW cartridge heaters, contained six wall thermocouples each. 

In measuring boiling heat transfer coefficients, great care must be 
exercised with the cartridge heater and temperature measuring 
instrumentation to ensure good accuracy. The instrumented test tubes were 
fabricated using the same method as that used by Hahne and Muller [Ref. 
17] and Wanniarachchi et al . [Ref. 14], The exact procedure can be found 
in Eraydin [Ref. 28]. Figure 7 is a cross-sectional sketch of an 
instrumented tube, showing the construction details and the location of 
the wall thermocouples. The thermocouples were embedded in the wall at 
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different circumferential and longitudinal positions along the heated 
section of the tube shown in Figure 7. 

The five instrumented tubes were located along the centerline of the 
tube bundle, forming a vertical in-line column. All the instrumented and 
active tubes were Turbo-B tubes made by Wolverine Tube Co. (see Section 
E) . These tubes were nominally 15.9 mm in outside diameter and were 
arranged in an equilateral triangular pitch (ie. centerline-to-centerline 
spacing) of 19.1 mm, giving a pitch-to-diameter ratio of 1.35. The 
thickness of the Turbo-B enhancement was 0.85 mm giving a diameter to the 
base of the enhancement of 14.2 mm. 

The bundle also contained a number of unheated dummy smooth tubes 
(marked "D") that were used to guide the two phase mixture through the 
bundle. The dummy tubes were made from commercially available 15.9 mm 0D 
smooth copper tubing. Two vertical baffle plates made of aluminum were 
used on either side of the bundle to restrict circulation into and out of 
the bundle at the sides. A dummy rack (Figure 8) consisting of 12 solid 
rods made of aluminum (15.9 mm 0D and spaced 19.1 mm from centerline-to- 
centerline) was placed below the tube bundle. This rack had a triangular 
pitch arrangement with vapor retaining plates on the sides and was 
designed for two purposes: to collect all rising two phase flow generated 
by the simulation heaters and direct it into the test bundle and to 
simulate vapor passing through a large bundle before reaching the 
instrumented tubes. A small open space (approximately 5 mm in height) was 
left between the bundle and dummy tube rack. This space allowed some 
refrigerant from below to enter the bundle and replace the vapor being 
generated in the bundle. However, there was also a space below the dummy 
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rack that allowed the majority of the circulation to occur. Thus, 
liquid/vapor circulation was vertically upward over the five instrumented 
test tubes with no net horizontal component. Most of the liquid-vapor 
mixture after passing through the bundle was separated when it reached the 
pool surface. However, due to the strong circulation patterns set up 
within the liquid pool, some vapor bubbles remained trapped in the liquid 
and circulated around the pool. 

The condenser included four test tubes (each 1.219 m in length and 
15.9 mm OD) in a vertical in-line column and five auxiliary copper coils. 
For the boiling experiments, these tubes were used to regulate the 
pressure and temperature in the evaporator. The condenser was designed to 
permit independent condensation studies of small in-line tube bundles, 
using the evaporator as a source of vapor. Details of the condenser can 
be found in Mazzone [Ref. 29] . 

D. DATA ACQUISITION SYSTEM/ INSTRUMENTATION 

As described by Akcasayar [Ref. 25] , a Hewlett Packard HP-3497A Data 
Acquisition System, HP-9125 computer and HP-701 printer were used for data 
acquisition, data reduction and data printing respectively. Although an 
HP-9826 computer and HP-7470A plotter can be used for final graph 
printing, a Macintosh computer (using Criketgraph) was utilized. As 
described by Anderson [Ref. 13], type-T copper-constantan thermocouples 
were used to make temperature measurements on the HP-3497A using a relay 
multiplexer assembly equipped with thermocouple compensation. A 20- 
channel relay multiplexer card was used to measure the voltage output from 
both the voltage and current sensors. The voltage measurements were taken 
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from separate sensors that measured the voltage going to the tube bundle, 
the simulation heaters and the auxiliary heaters. The total current 
going through the auxiliary and simulation heaters was measured using an 
American Aerospace Control (ACC) current sensor. The current to each 
instrumented tube heater was measured using five identical current 
sensors. The voltage supplied to the other active tubes was also 
measured, but the current of each active tube was not. Instead, the total 
current for a pair of active tubes was measured, and this was sufficient 
since these tubes each had the same power output (1000 W) as the 
instrumented test tube heaters and there was no apparent reason to monitor 
each active tube individually. Computer channel assignments for data 
acquisition and array assignments are given in Table 2. 

E. GEOMETRY OF TURBO -B TUBE 

The Turbo-B tube, manufactured by Wolverine Tube Inc., contained an 
enhanced surface geometry. The exterior boiling enhancement is made by 
raising low integral fins, cutting diagonally across these fins, and then 
rolling the fins to compress them to form mushroom-like pedestals [Ref. 
30]. This process forms numerous re-entrant passageways. Figure 9 shows 
the surface of the tube at 25 times its actual size. The tube is 
currently available in copper, cupro-nickel , and low carbon steel. 

The relative dimensions of the tube used in this study are as follows: 
Tube material - Copper 
Nominal Outside diameter — 15.9 mm 
Enhanced surface length — 203.2 mm 
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Thickness of Enhancement - 0.85 mm 
Diameter to Base of Enhancement * 14.2 mm 
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Figure 1. Schematic View of the Experimental Apparatus 
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Figure 2. Evaporator/Condenser Schematic 
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Figure 3. Front View of Evaporator 
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Figure 4. Side View of Evaporator 
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Figure 5. Photograph of 208 V, 75 A, Variable Transformers 



21 



X 




Figure 6. Sectional View of Evaporator Showing Tube Bundle 
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DUMMY TUBE RACK 





Figure 7. Thermocouple Locations on an Instrumented Boiling Tube and 
Tube Section View 
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Figure 8. Photograph of Dummy Rack 
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Figure 9. Close-up View of Turbo-B Tube Surface (25 X) 
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Table 1. EVAPORATOR HEATERS 



Heater Type 


Number 


Power Rating per Heater 


Instrumented Tube Heaters 


5 


1000W 


Active Tube Heaters 


12 


1000W 


Auxiliary Heaters 


4 


4000W 


Simulation Heaters 


5 


4000W 



Table 2. COMPUTER/DATA ACQUISITION 
ASSIGNMENT 



Amperage Sensor Description 


Channel 


Array 


Tube 1 


30 


Amp(O) 


Tube 2 


31 


Amp(1 ) 


Tube 3 


32 


Amp(2) 


Tube 4 


33 


Amp(3) 


Tube 5 


34 


Amp(4) 


Active Heater Group 1 


35 


Amp(5) 


Active Heater Group 2 


36 


Amp(6) 


Active Heater Group 3 


37 


Amp(7) 


Active Heater Group 4 


38 


Amp(8) 


Active Heater Group 5 


39 


Amp(9) 


Auxiliary Heaters 


25 


Amp(10) 


Simulation Heaters 


26 


Amp(1 1) 



Voltage Sensor Description 


Channel 


Array 


Instrumented/Active 


27 


Volt(O) 


Simulation Heaters 


28 


Volt(1) 


Auxiliary Heaters 


29 


Volt(2) 
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Table 2. COMPUTER/DATA ACQUISITION ASSIGNMENT (CONT . ) (cont . ) 



Thermocouple Description 


Channel 


Array in code 


Vapor 1*Top of Condenser 


00 


T(0) 


Vapor 2*Top of Condenser 


01 


T(1) 


Vapor 3-Top of Evaporator 


02 


T(2) 


Liquid 1-Top of bundle 


03 


T (3) 


Liquid 2-Top of bundle 


04 


T (4) 


Liquid 3-Bottom of bundle 


05 


T (5) 


Tube 1,No. 1 


40 


T (6) 


Tube 1,No. 2 


41 


T (7) 


Tube 1,No. 3 


42 


T(8) 


Tube 1,No. 4 


43 


T(9) 


Tube 1,No. 5 


44 


T(10) 


Tube 1,No. 6 


45 


TO 1) 


Tube 2, No. 1 


46 


T(12) 


Tube 2,No. 2 


47 


T( 1 3) 


Tube 2, No. 3 


48 


T ( 1 4) 


Tube 2,No. 4 


49 


T ( 1 5) 


Tube 2, No. 5 


50 


T( 1 6) 


Tube 2,No. 6 


51 


T(1 7) 


Tube 3, No. 1 


52 


T( 1 8) 


Tube 3,No. 2 


53 


T( 1 9) 


Tube 3,No. 3 


54 


T (20) 


Tube 3, No. 4 


55 


T(21) 


Tube 3, No. 5 


56 


T(22) 


Tube 3, No. 6 


57 


T(23) 


Tube 4, No. 1 


58 


T(24) 


Tube 4,No. 2 


59 


T (25) 


Tube 4, No. 3 


60 


T(26) 


Tube 4, No. 4 


61 


T(27) 


Tube 4, No. 5 


62 


T(28) 


Tube 4, No. 6 


63 


T (29) 


Tube 5, No. 1 


64 


T (30) 


Tube 5, No. 2 


65 


T (3 1 ) 


Tube 5, No. 3 


66 


T (32) 


Tube 5, No. 4 


67 


T(33) 


Tube 5, No. 5 


68 


T (34) 


Tube 5, No. 6 


69 


T (35) 
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IV. EXPERIMENTAL PROCEDURES 



A. REMOVAL OF THE TUBE BUNDLE AND BUNDLE DISASSEMBLY 

Before starting the removal of the tube bundle from the evaporator, 
the front glass viewing windows were carefully removed. Next, all 
thermocouple wires and tube heater electrical connections were 
disconnected. After this was completed, the nuts securing the backing 
plate and support block were removed, and the tube bundle was taken out 
from the back of the evaporator. 

When the bundle needed to be disassembled, it was ensured that there 
was a clean working surface. The first task was to remove the plexiglas 
plate attached to one end of the aluminum baffle plates (ie. the front of 
the bundle assembly) by four screws. The ten screws on the side of each 
aluminum baffle plate were then removed (these were attached to the dummy 
tubes down each side of the bundle) . The aluminum plates were then pulled 
off the bundle. The four corner dummy tubes (two top and two bottom) 
remained attached to the tube bundle support block as they were 
countersunk into the block. The six outer smooth tubes (three per side) 
could be easily pulled from the aluminum plates as they were attached only 
by the screws already removed. The other ten smooth tubes were then 
unscrewed from the tube bundle support block as seen in Figure 10. The 
smooth tubes were engraved to ensure proper identification during 
reassembly. With these tubes and aluminum baffle plates removed, only the 
instrumented and active enhanced heater tubes remained. These tubes were 
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removed from the support block by loosening the outer O-ring compression 
plate, disconnecting the active heater tube wired pairs, and pulling the 
tubes from the block. Reassembly of the tube bundle is done by reversal 
of this procedure. 

B. SYSTEM CLEAN-UP 

If the system had been previously operated with refrigeration oil (or 
contaminated from some other source) it had to be thoroughly cleaned. To 
accomplish this, the entire apparatus had to be taken apart and cleaned in 
the following manner. 

After removal of the refrigerant (R-113 by directly draining into 5 
gallon drums via drain valve R-5 (see Figure 2 in Chapter 3) at the bottom 
of the evaporator or R-114 by boiling off into the storage tank by opening 
R-l and R-7) and with the system at atmospheric pressure, all electrical 
connections to the bundle were disconnected and the front viewing glass 
windows were removed. An electric fan was used for safety to ensure 
proper ventilation. The tube bundle was then removed and disassembled as 
described in section A; the dummy tube rack was also removed. 

Having removed the tubes from the tube bundle, they were individually 
washed with warm water, rinsed and then wiped down with acetone. The 
smooth tubes were cleaned with Copper Brite (a commercial copper cleaning 
product) to remove any oxidation. They were also wiped down with warm 
water and then with acetone. The same procedure was followed for the 
Turbo-B tubes except they were not cleaned with Copper Brite for fear of 
clogging the channels. During the cleaning process, a soft bristled 
toothbrush was used to ensure the enhanced surface was cleaned properly, 
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exercising care not to interfere with the tube surface. The evaporator 
shell was cleaned in a similar manner, using warm water and acetone. 

C. INSTALLATION OF THE TUBE BUNDLE 

Once the tube bundle had been cleaned and reassembled (see section A) , 
and before tightening the backing plate nuts, the whole assembly was 
carefully guided back into the evaporator section, ensuring the plexiglas 
viewing cover of the tube bundle was not damaged. After the bundle was in 
position, it was ensured that the dummy tube rack was properly positioned 
below the bundle and that the vapor thermocouple positions were still 1.75 
cm above the bundle. Then, all the nuts were tightened equally on opposite 
sides to give equal compression on the gasket. To replace the front 
window, very small, equal torques (using a torque wrench) were applied 
circumferentially to each nut on the outer ring support in turn. After 
the window was in place, each tube (which extended through the outer 0- 
ring compression plate) was lightly tapped forward so as to touch the 
front-viewing window. The backing plate was then tightened and the 
individual tube 0-rings compressed, providing a good seal for the system. 
The compression plate had grooves for the tube 0-rings to sit in to help 
with proper alignment and ensure a good seal. 

D. SYSTEM LEAKAGE TEST 

After the system was isolated from the atmosphere and system integrity 
was restored, a Seargent Welch 10 SCFM vacuum pump was connected to the 
apparatus (via valves R-l and R-8) and the pressure taken down to 25 inHg 
vacuum. Valves R-l and R-8 were then secured and the system was left 
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untouched for at least 24 hours to see if there was any air leakage in. 
If there was significant leakage (>1 inHg over 24 hours), then the vacuum 
was broken by cracking open valve R-2 slowly (this ensured that no 
moisture entered the system) . The system was pressurized (with air) to 15 
psig through valve R-2. Large leaks could then be detected by simply 
listening to the air issuing from the system; small leaks were detected by 
spraying a soapy water solution to all surfaces where leaks were most 
likely to occur (front viewing glass gaskets, backing plate gasket, all 
fittings /valves coming off the condenser/evaporator, O-ring seals of the 
bundle tubes etc). Extreme care must be taken to ensure no moisture 
enters the inside of the heated tubes where the heater wires protrude. 
After all leaks were detected and corrected, the system was again 
subjected to a vacuum for another 24 hour period. If the vacuum held, 
then the system was ready to receive refrigerant. If not, the above leak 
correction test was repeated. 

E . REFRIGERANT 
1. Fill 

a. From System Storage Tank 

A refrigerant storage tank was used to store R-114 during 
modification/repairs to the system. The storage tank prevented discharge 
of the R-114 into the atmosphere and made the experimentation less costly. 
To fill the evaporator with R-114 from the storage tank, the ethylene 
glycol/water coolant temperature was first reduced to —15 °C. The system 
pressure was then maintained below the storage tank pressure (vapor 
pressure of R-114 at 20 °C is approximately 15 psig) by circulation of the 



31 



coolant through the condenser test tubes and auxiliary coils. Valves R-6 
and R-4 were then opened to draw the R-114 from the storage tank to the 
evaporator. The amount of refrigerant that was transferred was controlled 
by throttling valve R-6 to obtain the desired level in the evaporator. If 
required, additional R-114 could be transferred from a 68 kg storage 
cylinder to the system using valve R-2 (see section l.b). 
b. From Refrigerant Storage Cylinder 

To fill the apparatus from the 68 kg storage cylinder, the 
system pressure was reduced in the same way as above. A hose assembly 
containing a Drierite gas purifier was connected between the storage 
cylinder and valve R-2. A gas purifier was used not only to remove all 
impurities, but also to remove any water from the refrigerant. Once in 
place, both the storage cylinder valve and R-2 were opened until the 
desired refrigerant level was reached in the evaporator. 

2. Removal to the Storage Tank 

For tube replacement, system maintenance or system clean up, the 
R-114 was transferred to the storage tank. The ethylene glycol/water 
coolant temperature flowing through the storage tank was cooled to —15 °C; 
valves R-7 and R-8 were opened and the vacuum pump was turned on to put 
the storage tank under vacuum. Once the storage tank was under a 20 inHg 
vacuum, valve R-8 was shut and valve R-l was opened. Next R-114 was slowly 
boiled off to the storage tank by using the tube bundle, simulation and 
auxiliary heaters at a heat flux of 600 kW/m 2 (slow boiling is important 
to ensure minimum transfer of oil). As the refrigerant level decreased, 
individual heaters were turned off to ensure none were uncovered. Once 
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below the level of the heaters, the final few cm of R-114 was boiled off 
using heat from the atmosphere. Once all of the R-114 was transferred, 
valves R-l and R-7 were shut. 

F. OPERATION 

1* System Startup, Securing and Emergency Procedures 

See Appendix D 

2. Normal Operation 

The evaporator was filled with R-114 to a level of approximately 
10 cm above the top tubes in the bundle. Prior to operating the system, 
the 8 ton refrigeration unit was run for approximately an hour to reduce 
the ethylene glycol/water coolant in the sump to a temperature of —15 °C . 
The pressure in the evaporator/condenser was usually 12 to 15 psig if the 
system had been secured overnight. As the sump was brought to 
temperature, the data acquisition system and computer were turned on. 
This allowed the temperature in the system to be monitored during cool- 
down to saturation conditions. With this and pump number one running, 
one auxiliary condenser coil and the four condenser test tubes were used 
to bring the pool down to a subcooled condition (for R-114, approximately 
1 °C on all three pool thermocouples) . Subcooling of the refrigerant was 
done to ensure the pool had an evenly distributed temperature prior to 
starting a run. After reaching this subcooled condition, all coolant 
supply to the condenser was secured. The pool was then allowed to 'heat 
up' by conduction from the surroundings. Once a saturation temperature of 
2.2 °C was reached, the instrumented tube(s) (and simulation heaters for 
test 7) was/were switched on and set to the desired heat flux value. 
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This lengthy procedure was done to prevent the tubes from prematurely 
nucleating. The heat flux of the instrumented tubes was then slowly 
increased at desired intervals by adjusting the rheostat. For increasing 
heat flux, the data was taken with very small heat flux increments (every 
1000 kW/m 2 ) , waiting at least 5 minutes to attain steady state conditions 
at each heat flux. At all regions of the boiling curve (and especially 
near the onset of nucleate boiling) , two readings were taken at each heat 
flux to ensure accuracy. The bundle was continuously visually monitored 
through the observation windows. Figure 11 shows the tube bundle 
arrangements used during the experimentation. Test one was with only one 
instrumented tube turned on at any position within the bundle. Test two 
was with instrumented tubes one and two active in the bundle. Test three 
was with instrumented tubes one, two and three active in the bundle. Test 
four was with instrumented tubes one, two, three, and four active in the 
bundle. Test five was with all five instrumented tubes active in the 
bundle. Test six was all five instrumented tubes plus all five pairs of 
active enhanced heater tubes active in the bundle. Test seven was the 
same as test six with the addition of all five simulation heaters active. 
For each data set, the five simulation heaters had the same heat flux as 
the tubes within the bundle. 

G. OIL ADDITION 

During the bundle experiments, successive amounts of York-C oil were 
added into the evaporator. Since the weight of the refrigerant in the 
evaporator was 60.3 kg, the amount of oil corresponding to 1% by weight 
was measured as 670 ml, 2% 1340 ml etc. The oil was syphoned into the 
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evaporator via a funnel/hose connection through valve R-3 by reducing the 
pressure in the evaporator to less than 15 inHg vacuum. Ensuring that no 
air entered the system, valve R-3 was promptly shut when the desired 
amount of oil had been added. 

H. DATA REDUCTION PROCEDURES 

The data reduction program "DRP4RH" was used during the experiments 
for processing the data collected ( see Appendix E for listing) . The 
program was written in HP Basic 3.01 and run on an HP-9000 series 
computer. The characteristics and capabilities of this software are 
similar to those provided by Anderson [Ref. 13] . The following 
modifications were made: 

I. Correction for pool height by Chilman [Ref. 3] 

2. Installation of new thermocouple at the bottom of the liquid pool 
(bundle inlet temperature) by Chilman [Ref. 3] 

3. The ability to obtain data from one instrumented tube at any position 
in the bundle . 

4. The natural convection correlation of Churchill and Chu [Ref. 31] for 
a single horizontal cylinder in an 'infinite' liquid pool was added for 
comparison with experimental data. 
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Figure 10. Photograph of Tube Bundle Support Block 
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Figure 11. Tube Bundle Arrangements 



used During Experimentation 



Test Four Test Five Test Six Test Seven 





V. RESULTS AND DISCUSSION 



A. INTRODUCTION 

The results are presented in four sections with subs-sections as 
appropriate. The first section discusses the preliminary experiments 
which led to modifying the experimental start-up procedure to include 
subcooling. The second section discusses the natural convection effects, 
nucleate pool boiling phenomena and hysteresis effects within the tube 
bundle in pure R-114. The third section discusses similar phenomena for 
R-114/oil mixtures and their effects on the above. The fourth section 
shows comparisons of data taken during this thesis with previously 
obtained data at the Naval Postgraduate School. 

A list of data files taken during this investigation may be found in 
Appendix A. All data files used in this thesis use the following filename 
sequence. Each file is composed of five sets of alpha-numeric characters 
used to describe the experiment. 



First set (2 char.) Tube Type 



TB (Turbo-B) 



Second set (1 char.) Heat Flux 



I (Increasing) 
D (Decreasing) 



Third set (2 char.) Oil Percent 



00 ( 0 %) 
01 ( 1 %) 
02 (2%) 
03 (3%) 
06 (6%) 
10 ( 10 %) 



Fourth set (2 char.) Test Type 



01 (test 1) 

02 (test 2) 

03 (test 3) 

04 (test 4) 
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05 (test 5) 

06 (test 6) 

07 (test 7) 



Fifth Set (1 char.) Additional tests A-Z (If conducted) 

To give an example, the filename TBI0107 means "Turbo-B tube bundle, 
increasing heat flux, 1% R-114/oil mixture and test number 7". If more 
detail is desired about a specific data set, see Appendix A. All plot 
filenames are similar to data filenames except they start with the letter 
M P M . The test numbers are shown in Figure 11 in Chapter IV. 

All graphs are plotted showing heat flux (W/m 2 ) along the ordinate (y 
axis) as a function of wall superheat (K) along the abscissa (x axis). 
The wall superheat is defined as the difference between the 'corrected' 
average tube wall temperature (ie. having accounted for depth of 
thermocouple burial) and the local liquid saturation temperature 
(corrected for hydrostatic head within the bundle) . The heat flux was 
corrected to account for the heat lost through the unheated tube ends. 
The heat flux was varied from 600 to 100,000 W/m 2 for increasing heat flux. 
To ensure greater detail at the point of incipience, the heat flux was 
increased in small steps; these settings varied from test to test. The 
heat flux values for decreasing experiments were taken at prescribed 
settings for easy comparison with past experiments and future reference. 
These heat flux settings were 1x10 s , 7.5x10*, 5x10*, 3x10*, 2x10*, 1.5x10*, 
1x10*, 7xl0 3 ,4xl0 3 , 2xl0 3 , lxlO 3 W/m 2 . Approximately 30-40 data points were 
taken for each increasing heat flux run and 20-25 points were taken for 
each decreasing heat flux run. 
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B . PRELIMINARY EXPERIMENTS 



After the pure R-113 was removed following Chilman's [Ref. 3] 
experiments and the apparatus and system cleaned, pure R-114 was added 
from the storage cylinders. Five tests were conducted using this R-114 
and these are shown in Figure 12. The first test (TBI0001A) was test one 
with the top tube activated. The procedure described by Eraydin [Ref. 28] 
was followed, but the plot of test TBI0001A and Eraydin' s data for test 
one (also shown in Figure 12) produced significantly different results in 
the natural convection (NC) region. The data of Eraydin show a greater 
heat transfer coefficient (lower wall superheat) than test TBI0001A which 
show results closer to the Churchill/Chu [Ref. 31] correlation (C/C) for 
natural convection. The only difference in the apparatus between test 
TBI0001A and Eyradin's experiments was the addition of a third 
thermocouple at the bottom of the pool. For test TBI0001A, all three pool 
temperature readings were within +/- 0.1 °C prior to recording data. For 
the data of Eraydin, only the temperature at the top of the pool could be 
checked. 

Test TBI0001B was thought to be a repeat of test TBI0001A, but for the 
bottom tube in the bundle (tube 5 only) . However, upon observation and 
investigation, it was found that two tubes (tube one and tube five) were 
activated due to the way the program DRP4 was set up. Hence the data 
presented is for tube one (the top tube) with tube five (the bottom tube) 
activated as well. The program was then modified to obtain data for a 
single tube (test one) at any position within the bundle. 

Test TBI0001C was conducted using the bottom tube (tube five) as a 
single tube. Again following Eraydin' s [Ref. 28] procedure, partial 
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nucleation was observed immediately, explaining why the data lie well to 
the left of the Churchill/Chu correlation (ie. a higher heat transfer 
coefficient) . It can be seen that test TBI0001C is in good agreement with 
the single tube of Eraydin (top tube). The reason for this is probably 
due to the start-up conditions, which were not carefully monitored for any 
of these tests. 

Test TBI0001D was conducted using only tube three and also followed 
the procedure of Eraydin. This displayed similar behavior as test 
TBI0001A (tube one) except that nucleation was delayed, ie . occurred at a 
higher wall superheat. 

All experiments thus far were conducted following Eraydin' s procedure 
for pure R-114. It was next decided to vary his procedure slightly. For 
test TBI0001E, the pool was first subcooled slightly to 1°C . This ensured 
that any nucleation sites were deactivated and that the whole pool was at 
an even temperature. The pool was then brought up to the required 
saturation temperature. However, partial nucleation was still observed 
during the run. The plot of TBI0001E is similar to test TBI0001C and 
Eraydin' s test one. 

With all of these confusing results, it was decided to empty, clean 
and recharge the evaporator with fresh R-114. Upon boiling off the 'old' 
R-114, a small quantity of oil contamination was found in the bottom of 
the evaporator. This could have been either vacuum pump oil (which may 
have leaked into the apparatus) or refrigerant oil that entered the system 
with the R-114 charge (which should be minimal) . A third possibility 
could be oil from previous refrigerant/oil mixture experiments. However, 
this is not likely since the system had been completely stripped and 
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cleaned twice since the last mixture tests had been conducted. Samples of 
each type of oil (vacuum pump oil, miscible refrigerant oil (York-C) and 
a sample of the contamination) were sent away for analysis and the results 
were still forthcoming at the time of writing this thesis. From the color 
of the contamination, it would appear to have originated from the vacuum 
pump. This would seem to be more logical since Chilman had experienced 
problems with the vacuum pump. 

The system, including the bundle, were cleaned thoroughly and fresh 
R-114 was added. However, a gas purifier was utilized to ensure that only 
pure, clean R-114 (with no moisture) was added (see Chapter IV for further 
details). Test TBI0001F then was conducted following Eraydin's [Ref. 28] 
procedure with no subcooling and 'ignoring' the bottom pool thermocouple 
value. The data showed good agreement with Eraydin in the natural 
convection region, but a lower heat-transfer coefficient in the boiling 
region. A possible explanation for this lower heat transfer is the fact 
that each individual tube was fully cleaned prior to adding the new R-114 
and the surface characteristics may have been modified in some way. 

Test TBI0001G was conducted using the same procedure except the pool 
was initially subcooled down to 1°C for 30 minutes. This ensured the pool 
had an even temperature distribution throughout. The data were now much 
closer to the Churchill/Chu [Ref. 31] correlation. However, some 
premature nucleation was still observed. Test TBI0001J repeated the above 
test with the pool temperature subcooled to 1°C for at least one hour to 
further deactivate any remaining nucleation sites within the bundle prior 
to starting experiments. The data then agreed with the Churchill/Chu 
correlation as seen in Figure 12. It became apparent that premature 
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nucleation could affect the natural convection data significantly. 
Therefore, for all subsequent tests, this same procedure was adopted with 
strict observation of the bundle to ensure no premature nucleation 
occurred. 

C. PURE R-114 TURBO -B TUBE BUNDLE EXPERIMENTS 
1. Test One for Different Tube Positions 

The first set of experiments conducted were performance tests in 
pure R-114. Figure 13 shows increasing heat flux for a single tube within 
the bundle at different positions (positions 1, 3, and 5) while Figure 14 
shows corresponding data (including typical uncertainties) for decreasing 
heat flux. All three tube positions agree closely with the Churchill/Chu 
[Ref. 31] correlation in the natural convection region. The difference 
between position 1 and position 3 and 5 may be that position 1 is affected 
by the fact that the flow is 'free' to expand after leaving the bundle. 
This difference may also be due to wall temperature uncertainty due to 
differences in the fabrication process (see uncertainty analysis Appendix 
C) ; however in the natural convection region, this uncertainty is low due 
to relatively high values of wall superheat. Figure 13 also shows that 
tube position within the bundle may influence the point of incipience. 
Bergles and Rohsenow [Ref. 32] have studied the incipient point in more 
detail. They concluded that nucleation was controlled by a nucleation 
parameter, N, given by 

(*)(T s )(v fg ) 

N - 

(Pv)(hfg) 
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which creates an incipient boiling superheat given by 
(T w - T s ) - 2(N)/r. 

In the above expression, r is the local bubble radius. Calculation of the 
nucleation parameter for pure R-114 showed that as saturation pressure 
increases, the nucleation parameter decreases. Assuming that the radius 
of curvature of the forming bubbles is constant (which is reasonable for 
a Turbo-B surface which has large, regularly spaced cavities), then (T g — 
T s ) also decreases and nucleation may be expected to occur earlier. This 
certainly seems to be verified in Figure 13 where the incipient point 
occurs earlier (lower wall superheat) for a lower tube ie. where there is 
an increase in the local saturation pressure. You et al . [Ref. 33] also 
showed a decrease (approx. 30%) in the average incipient superheat as 
pressure was increased from 1 to 1.5 bar for pool boiling of FC-72 on a 
single tube, offering some other experimental verification for this 
conclusion . 

Chilman [Ref. 3] conducted test one for R-113 using the top tube 
only and varied the local saturation pressure by varying the pool height 
in the evaporator. He found that the point of incipience was delayed when 
increasing the liquid pool height (ie. the hydrostatic pressure head). 
Boundary layer effects due to different liquid circulation patterns may 
have caused this delay in nucleation and more research is certainly needed 
to fully understand the influence of pressure on the point of incipience. 
One experiment that could be conducted would be to vary the pool height, 
but keep the local pressure at each tube constant by simultaneously 
varying the vapor pressure above the pool. 
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Once nucleation occurs, Figure 13 shows that the single tube 
experiments merge onto a single boiling curve. In this region, there 
appears to be no effect of hydrostatic pressure head. Figure 14 shows the 
corresponding decreasing heat flux data for a single tube at the same 
three positions within the bundle. It shows no significant influence of 
tube position in the bundle for decreasing heat flux. Note that at low 
heat flux, the experimental uncertainty in heat flux and wall superheat is 
the largest (see Appendix C) . 

2. Test Two to Test Seven 

Figures 15 to 20 show test two to test seven for increasing heat 
flux with pure R-114. Also shown for comparison in each figure is the 
Churchill/Chu [Ref. 31] correlation although this is only truly valid for 
a single tube in an infinite pool. Figure 15 shows good agreement with 
the Churchill/Chu correlation in the natural convection (NC) region and 
shows no effect of the lower tube on the upper tube. The incipient point 
occurs approximately at the same wall superheat for both tubes; once 
boiling the lower tube has the higher heat-transfer coefficient. This is 
contrary to the results obtained by Chilman [Ref. 3] and Anderson [Ref. 
13] with pure R-113, where the higher tube had the better heat transfer. 
The reason for this difference is not known, but may be due to the 
explosive nature of incipience for R-114 compared to the more gradual 
partial incipience for R-113. For R-114 experiments, the pool was 
subcooled by 1 °C for only about hour while for the R-113 experiments, the 
pool was left in a subcooled state since the previous experiment (for R— 
113 experiments, the pool is 'heated up' to saturation conditions). This 
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difference in subcooling may significantly affect the nature of the 
observed incipience. Further research should be conducted at the 
incipient point to address some of the questions. Figure 16 for three 
tubes activated shows similar behavior as Figure 15 (ie. during nucleate 
boiling, the lowest tube has the highest heat-transfer coefficient). 
Also, Figure 16 shows that the lowest tube seems to nucleate last. 
Figures 17 and 18 for four and five tubes activated show similar trends in 
both the NC and boiling regions. It appears that the tubes nucleate in 
order down the bundle (ie. the top tube nucleates at the lowest wall 
superheat and the bottom tube nucleates at the highest wall superheat) . 
For tests six and seven (Figures 19 and 20 respectively) , the maximum 
controllable heat flux was less than tests one through five due to the use 
of smooth tubes in the condenser limiting the condensate rate (and hence 
pressure) in of the vapor space (if enhanced tubes had been used, this 
could have been increased) . Figure 19 (test six) is consistent with the 
above trends. Furthermore, the effect of activating the whole bundle 
seems to cause the lowest tube to nucleate at a lower heat flux and wall 
superheat. It also appears that there may be some influence of tube 
position in the boiling region; however, this is probably due to 
inaccuracy in the wall temperature measurements (see Appendix C) . When 
the simulation heaters are also activated (test seven, Figure 20) the 
trends in the NC region are similar (ie. no effect of lower tubes on upper 
tubes). Full nucleation of the bundle however, seems to occur earlier. 

Figure 21 shows the data from tube one for all seven tests with 
increasing heat flux. This figure is of more fundamental interest as it 
shows the same tube under different bundle conditions. It therefore gives 
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a better direct comparison of the effect of heated lower tubes as any 
uncertainty in the tube wall measurements (the largest error in the 
experimental data) is the same for each test (ie. any effects seen in the 
data are bundle effects). In the NC region, tube one alone (test one) is 
somewhat different. This may be due to 'expansion' of the flow as it 
leaves the top of the bundle (ie. where the velocity of the flow has 
slowed down) or due to convective effects by the addition of another tube. 
For test two, Figure 21 shows an effect of the lower tube on tube one 
performance in the NC region. For test three and all subsequent tests no 
further improvement is seen. It appears, therefore that in the natural 
convection region, an upper tube is affected by a lower tube directly 
below; however, when additional lower tubes are heated, there is no 
further increase in performance of the top tube. There is also no effect 
on the incipient point (apart from test one mentioned above) . In the high 
heat flux boiling region, there is also little enhancement due to the 
lower tubes. This is to be expected in an enhanced tube bundle, where the 
total heat transfer at high heat fluxes is primarily due to nucleation 
from the tube surface itself, rather than from convection around the 
surface from the tubes below. 

Figures 22 to 27 show tests two to seven for decreasing heat flux 
in pure R-114. When comparing tube one to tube two, Figure 22 shows a 
increase in heat transfer performance of tube one by tube two in the 
boiling region while Figure 15 (increasing heat flux) showed the opposite 
effect. The most probable reason for this crossover is that these two 
experiments (TBI0002) and (TBD0002) were conducted on different days and 
startup procedures were slightly different. For TBI0002 (increasing heat 
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flux), the test was conducted as outlined in Chapter IV section F (ie. 
subcooled to 1 °C , gradually heated up with data taken over a period of 
approximately 4 hrs) . For TBD0002 (decreasing heat flux), the pool was 
subcooled to 1 °C , and then tube one and two were turned on to the highest 
heat flux (10 5 W/m 2 ) and allowed to heat up for 30 minutes. The total time 
boiling for TBI0002 at the highest heat flux was therefore less than for 
TBD0002. For all other experiments, increasing followed by decreasing 
runs were conducted on the same day approximately 15 minutes apart. More 
research should be conducted to investigate nucleation site 
activation/deactivation. 

At high heat fluxes for test three, Figure 23 shows no heat 
transfer performance improvement between the tubes. However, at low heat 
fluxes, there does appear to be an improvement on tube one and two from 
tube three. Figures 24 to 27 show similar trends at high and low heat 
fluxes (ie. the top tubes are further enhanced by lower tubes at low heat 
fluxes) . It should also be noted that the lowest tube in any specific 
decreasing heat flux test had the worst performance. 

Figure 28 compares tube one only for tests one to seven for 
decreasing heat flux. As stated above, there appears to be a definite 
tube enhancement at low heat fluxes with little or no enhancement at high 
heat fluxes. This is probably due to convective effects which tend to 
increase the heat transfer performance of the upper tubes due to the 
presence of lower tubes (ie. bubbles coming from the lower tubes impinge 
and slide over the upper tubes and increase the heat transfer) . At high 
heat fluxes, on the other hand, all the tubes are nucleating so vigorously 
that these ' sliding 1 bubbles have little or no noticeable affect on the 
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overall performance. This supports the hypothesis of Cornwell [Ref. 7] 
that total heat transfer in a bundle is due to a summation of convective 
and nucleation heat transfer phenomena. Similar trends were found by 
Anderson [Ref. 13] and Akcasayar [Ref. 25] for smooth and finned tube 
bundles respectively (using the same apparatus) and also by Arai et al. 
[Ref. 20] for a Thermoexcel-E tube bundle. However, Akcasayar [Ref.] did 
not find such an enhancement effect for a High Flux tube bundle indicating 
that at low heat fluxes, a porous coated surface already has a significant 
number of active nucleation sites such that impinging bubbles from below 
have little or no added effect. Turbo-B is more similar to a Thermoexcel- 
E surface and at low heat fluxes , these two types of surface obviously 
exhibit different nucleation characteristics to those of a porous coated 
tube . 



D. R-114/0IL MIXTURES TURBO-B TUBE BUNDLE EXPERIMENTS 
1. Tests with 1% and 2% oil 

Only four experiments were conducted with a 1% and 2% R-114/oil 
mixture. These were tests one and seven for both increasing and 
decreasing heat flux; no experiments were conducted for tests two through 
six. Figure 29 shows test one at 1% oil concentration for increasing and 
decreasing heat flux, clearly showing a hysteresis pattern. Compared with 
pure R-114 (Figures 13 and 14) , Figure 29 shows no apparent effect of oil 
on the heat transfer in either the NC or boiling regions. 

Figure 30 shows test seven for increasing heat flux (1% oil) . 
Compared with Figure 20 (pure R-114), there are similar trends (ie. no 
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effect in the NC or boiling regions) . The tubes again appear to be 
nucleating ' in order' down the bundle, as found with test seven in pure R— 
114. For decreasing heat flux, Figure 31 shows similar trends to test 
seven in pure R-114 (Figure 27) . Thus a 1% oil concentration appears to 
have little or no effect on bundle performance for both increasing and 
decreasing heat flux. 

Figure 32 shows test one for a 2% oil concentration for 
increasing and decreasing heat flux. Figures 33 and 34 show test seven 
for increasing and decreasing heat flux respectfully for 2% oil 
concentration. All three graphs (Figures 32 to 34) are similar to those 
for pure and 1% oil concentrations, showing that 2% oil also has little 
effect on overall bundle heat transfer performance. 

2. Tests with 3% oil 

Nine experiments were conducted with a 3% R-114/oil mixture. In 
addition to test one and seven (conducted for both increasing and 
decreasing heat flux as before), tests two through six were conducted for 
decreasing heat flux only. Figure 35 shows test one with 3% oil for both 
increasing and decreasing heat flux. Again, the figure clearly shows 
hysteresis effects between the increasing and decreasing experiments. As 
with previous oil percentages (Figures 13, 14, 29, and 32), it shows there 
is no apparent effect of oil in the NC region, similar to previous test 
one data for other oil percentages. However, at the highest heat flux 
(100 kW/m 2 ) there is an increase in the heat transfer of about 10%. This 
is similar to the increases found by Burkhardt and Hahne [Ref. 23] in a 
finned tube bundle and Arai et al. [Ref. 20]. Figure 36 shows test seven 
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for increasing heat flux. The Churchill/Chu [Ref. 31] correlation for 
pure refrigerant is plotted for comparison only. Agreement is good, 
demonstrating that in the NC region, in addition to there being no effect 
of tube position, there is also no apparent effect of oil concentration on 
the heat-transfer coefficient. As before, the tubes appear to be 
nucleating ' in order' (ie. top tube nucleates first with the bottom tube 
nucleating last) . 

Figures 37 to 42 show data from tests two to seven for decreasing 
heat flux only. All show no effect of lower tubes on upper tubes (within 
the bundle) in the boiling region at high heat fluxes (the small amount of 
scatter is probably due to inaccuracies in the wall temperature 
measurements). Each successive figure shows that the lowest tube has the 
lowest heat-transfer coefficient; this tube is then enhanced by the 
activation of tubes below it. Again it should be noted that the 
experimental uncertainty is larger at low heat fluxes. At all oil 
concentrations, tube five is seen to have the lowest heat transfer 
performance. According to Chilman [Ref. 3], tube five had the highest 
uncertainty in the wall temperature measurements and this might be the 
cause of this discrepancy. 

If one compares Figure 42 for 3% oil with Figure 28 for pure R- 
114, it can again be seen that there is a small increase in the bundle 
heat-transfer coefficient for the R-114/oil mixture at the highest heat 
fluxes. For all tests with oil added, significant foaming was observed at 
the pool surface and this may be the cause of this increase in heat 
transfer. Schlager et al . [ref. 21] in their review article point out 
that for certain conditions (typically low pressure and high heat flux), 
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the heat-transfer coefficient increases at low oil concentrations; they 
attributed this to foaming. Figure 43 compares test one to seven for tube 
one for decreasing heat flux. As before with pure R-114 (Figure 28), 
there appears to be a definite increase in performance of the upper tubes 
by lower tubes at the low heat fluxes due to convection effects, with 
little or no such increase at high heat fluxes. 

3. Tests with 6% oil 

The same nine tests as with 3% oil were conducted with a 6% R- 
114/oil mixture. Figure 44 shows test one for both increasing and 
decreasing heat flux. It clearly shows a hysteresis 'loop' between 
increasing and decreasing experiments. In comparison with other oil 
concentrations, the point of incipience occurs at a slightly lower heat 
flux. There also appears to be a small degradation in performance (10-15% 
compared with 3% oil concentration) at the highest heat flux (100 kW/m 2 ) 
due to the oil, but there is no apparent effect in the NC region. Figure 
45 shows test seven for increasing heat flux. As before, there is no 
apparent effect of oil on the heat transfer in the NC region and the tubes 
appear to be nucleating 'in order'. The point of incipience also seems 
unaffected by the presence of the oil. 

Figures 46 to 51 show tests two to seven for decreasing heat 
flux. At the highest heat fluxes, there is a similar small degradation 
in the heat transfer as found with test one (10-15%) when compared with a 
3% oil concentration (Figures 37 to 42) . When compared with pure R-114 
(Figures 22 to 27), there is neither enhancement nor degradation, 
indicating that any enhancement provided by 3% oil is offset by 6% oil. 
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At low heat fluxes, the data are not only very similar to that for pure 
refrigerant, but also to the other R-114/oil mixtures (ie. at low heat 
fluxes, there is no effect on heat transfer at any oil concentration). 
Figure 52 compares tests one to seven for tube one for decreasing heat 
flux. As before, there appears to be the same convective enhancement at 
low heat fluxes with no enhancement (due to the successive activation of 
lower tubes within the bundle) at high heat fluxes. 

4. Tests with 10% oil 

The same nine tests were repeated for an R— 114/oil mixture with 
10% oil. Figure 53 shows test one for both increasing and decreasing 
heat flux. Incipience occurred at a slightly heat higher flux than both 3% 
and 6% oil concentrations indicating that there appears to be no 
systematic increase or decrease in this point with increase in oil 
concentration. More importantly, there is a significant decrease in the 
heat transfer at the highest heat fluxes (20%) when compared with pure R— 
114. This is probably due to the re-entrant channels becoming 'clogged' 
with oil as the R-114/oil mixture is 'transported' to the surface at a 
high rate. Figure 54 shows test seven for increasing heat flux. As 
before, this shows that the NC region is unaffected by either oil 
concentration or lower tubes in the bundle. At the highest heat fluxes 
available (40 kW/m 2 ) there appears to be little decrease in the bundle 
performance (when compared to pure R-114) due to the oil. This indicates 
that at typical evaporator operating heat fluxes, the presence of oil does 
not significant effect the heat transfer enhancement process. At higher 
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heat fluxes, however, the effect of oil appears to be very significant as 
seen in Figure 53. 

Figures 55 to 60 show test two to test six for decreasing heat 
flux with 10% oil. At low heat fluxes, there seems to be no effect of the 
oil on the local heat transfer performance. However, at high heat fluxes, 
there is a significant decrease in performance. Interestingly, if one 
compares Figures 53, and Figures 55 to 58 at high heat fluxes, the lowest 
activated tube in the bundle is significantly degraded. The effect of 
activating a lower tube significantly enhances the heat transfer from the 
tube directly above and (to a lesser degree) the tubes even higher in the 
bundle. This may be due to the vigorous boiling action of lower tubes 
partly 'scouring' the oil rich layer which 'blankets' the upper tubes. 
This effect was also noticeable with the High Flux bundle (Akcasayar [Ref. 
25]). 

Figure 61 compares tube one for tests one to seven for decreasing 
heat flux. If one compares Figures 26 (0%), 43 (3%), 52 (6%) and 61 (10%) 
for tube one for all seven tests, it is clear that at low heat fluxes, the 
heat transfer coefficient is similar, regardless of oil concentration. 
Furthermore, convective effects are consistent and provide similar 
enhancements in heat transfer performance for all concentrations. At high 
heat fluxes at 0, 3, and 6% oil concentrations, there is little 
enhancement due to activation of lower tubes. However, at 10%, there does 
appear to be a small heat transfer enhancement due to activation of lower 
tubes. This was attributed above to increased 'scouring' of the oil from 
the vicinity of the Turbo-B surface by the increase in bubble activity as 
more tubes are activated within the bundle. However, for a practical 
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operating heat flux range between 15 and 30 kW/m 2 , there is no significant 
degradation in heat transfer performance for an oil concentration of up to 
10 %. 

E. COMPARISON OF R- 114/OIL MIXTURE EXPERIMENTS 

Figures 62 and 63 compare tube one from test one for increasing and 
decreasing heat flux for all oil concentrations. Figure 62 shows no 
effect of oil in the NC region, but some degradation in the boiling 
region. The correlation of Churchill/Chu is included for comparison. The 
incipient point appears relatively random, indicating no early or delayed 
nucleation caused by the presence of oil. Figure 63 shows similar 
degradation with a significant effect of the oil (20% decrease in the heat 
transfer from 0% to 10% oil) at the highest heat fluxes. 

Figures 64 and 65 compare tube one from test seven for increasing 
and decreasing heat flux for all oil concentrations. As compared to tube 
one test one (Figures 62 and 63) , the presence of oil has some degradation 
effect on the heat transfer performance (15%) in the NC region; this may 
be due to a change in the mixture properties which would tend to increase 
the wall superheat slightly as shown. As with Figure 62, Figure 64 shows 
that the incipient point appears relatively random. Figure 65 shows 
similar trends to Figure 63 (ie. no effect of oil at high heat fluxes), 
but shows convective effects at low heat fluxes. This was expected and 
previously reported (Figures 28, 43, 52, and 61). 

Figures 66 and 67 show the average bundle heat-transfer coefficient 
(ie. an average of all five instrumented tubes) as a function of heat flux 
for test seven for increasing and decreasing heat flux respectively at all 
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oil concentration. The data are from the same data set as that shown in 



Figures 64 and 65. Comparing Figures 64 and 66, degradation is seen in 
the NC region (15%) due to the change in mixture properties as mentioned 
above. However, at a practical operating heat flux range between 15 and 
30 kW/m 2 , the presence of up to 10% oil causes no degradation in bundle 
performance as seen. These trends over this heat flux range were similar 
to that found for a High Flux tube bundle (Akcasayar [Ref. 25]). 
Comparing Figures 65 and 67, similar trends (ie. no significant effect of 
oil) are found. However, due to the limit in maximum controllable heat 
flux for test seven, data at 'higher' heat fluxes ( up to 10 5 W/m 2 ) could 
not be obtained (as mentioned earlier). As shown in Figures 62 and 63, 
there may be a significant degradation in the heat transfer performance at 
these higher heat fluxes. 

F. COMPARISON WITH PREVIOUS NPS DATA 

Figure 68 shows a comparison between the present data for a Turbo-B 
bundle, the data of Anderson [Ref. 13] for a smooth tube bundle and the 
data of Akcasayar [Ref. 25] for both a 19 fpi and High Flux tube bundle in 
R-114. For clarity, only test one (tube one) for a decreasing heat flux 
in pure R-114 has been shown. Figure 68 shows that the Turbo-B tube has 
a significantly lower heat-transfer coefficient than the High Flux tube at 
all heat flux. This is surprising since Sugiyama [Ref. 34] showed that in 
the single tube apparatus, the Turbo-B tube was the best performer. The 
reason for this difference in behavior is not known. Also, all of the data 
appear to be parallel to each other; one may expect the enhanced tubes to 
have a different slope to a smooth tube due to the greater amount of 
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nucleation. At the highest heat fluxes, Turbo-B and the 19 fpi tube 
appear to have a similar heat transfer performance. However, the finned 
tube heat flux is based on the root diameter; if the actual finned area 
had been used, then the heat flux would be significantly lower. Figure 68 
shows that the heat transfer enhancement given by the Turbo-B tube when 
compared to a smooth tube is about three at high heat fluxes and increases 
to about five at low heat fluxes. 

Figures 69 and 70 compare the average bundle heat-transfer coefficient 
(ie. test seven) for a given oil percentage to that with no oil for all 
four tube bundles at heat fluxes of 15 and 30 kW/m 2 respectively. 15 and 
30 kW/m 2 were chosen as being representative of the lower and upper limits 
of heat fluxes used in practical Naval evaporators. At 15 kW/m 2 , Figure 
69 shows large enhancements for the smooth and finned tube bundles for all 
oil concentrations, especially at the lower oil concentrations. However, 
the Turbo-B and High Flux tube bundles show a degradation in the heat- 
transfer coefficient at all oil concentrations (approximately 5-10% at low 
oil concentrations dropping to nearly 25% at 10% oil for the High Flux 
bundle) . At the higher heat flux (30 kW/m 2 ) , Figure 70 shows similar 
trends as Figure 69. However, the High Flux bundle now exhibits a 40% 
decrease in the average bundle heat-transfer coefficient at 10% oil. 

Figures 71 and 72 compare the average bundle heat-transfer coefficient 
(ie. test seven) for each enhanced tube to that for the smooth tube bundle 
(tested by Anderson [Ref. 13]) for all oil concentrations at heat fluxes 
of 15 and 30 kW/m 2 respectively. At 15 kW/m 2 with pure R-114, Figure 71 
shows an enhancement factor of 3.7 for the Turbo-B tube bundle. This 
enhancement decreases slowly with increasing oil percentage to a factor of 
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about 2.5 at 10% oil. This agree very closely with the 19 fpi bundle. 
The High Flux bundle exhibits much larger enhancements , from over 6 at 0% 
oil to just over 3 at 10% oil. At 30 kW/m 2 with pure R-114, Figure 72 
shows an enhancement of 3.8 for the Turbo-B tube bundle, decreasing to 
about 3 at 10% oil. This again agrees closely with the 19 fpi bundle. It 
should be noted that the High Flux bundle enhancement has decreased to a 
similar value and gets worse than the other bundles as the heat flux is 
further increased. 
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F igure 12. Performance of Test One For Preliminary Experiments 
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Figure 13. Performance of Test One at Various Tube Positions for 
Increasing Heat Flux 
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Figure 14. Performance of Test One at Various Tube Positions for 
Decreasing Heat Flux 
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Figure 15. Performance of Tubes 1 and 2 for Increasing Heat Flux in Pure 
R-114 
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Figure 16. Performance of Tubes 1, 2, and 3 for Increasing Heat Flux in 
Pure R— 114 
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Figure 17. Performance of Tubes 1, 2, 3, and 4 for Increasing Heat Flux 
in Pure R-114 
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Figure 18. Performance of All Five Tubes for Increasing Heat Flux in 
Pure R-114 



65 



100000 



o 




O O O 

O O O 

O O <«- 

O <r- 

(31N/M) xmd 1V3H 

Figure 19. Performance of All Five Tubes with Active Pairs for 
Increasing Heat Flux in Pure R-114 
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Figure 20. Performance of the Bundle with Simulation Heaters for 
Increasing Heat Flux in Pure R-114 
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Figure 21. Comparison of Tests One to Seven for Tube 1 for Increasing 
Heat Flux in Pure R-114 
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Figure 22. Performance of Tubes 1 and 2 for Decreasing Heat Flux in Pure 
R-114 
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Figure 23. Performance of Tubes 1, 2, and 3 for Decreasing Heat Flux in 
Pure R— 114 
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Figure 24. Performance of Tubes 1, 2, 3, and 4 for Decreasing Heat Flux 
in Pure R-114 



71 



100000 



o 

o 




o 'r- 

T— 

(2IAI/M) xmd 1V3H 

Figure 25. Performance of All Five Tubes for Decreasing Heat Flux in 
Pure R-114 
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Figure 26. Performance of All Five Tubes with Active Pairs for 
Decreasing Heat Flux in Pure R-114 
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Figure 27. Performance of the Bundle with Simulation Heaters for 
Decreasing Heat Flux in Pure R— 114 
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Figure 28. Comparison of Tests One to Seven for Tube 1 for Decreasing 
Heat Flux in Pure R-114 
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Figure 30. Performance of the Bundle with Simulation Heaters for 
Increasing Heat Flux in R-114 with 1% Oil 
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Figure 31. Performance of the Bundle with Simulation Heaters for 
Decreasing Heat Flux in R-114 with 1% Oil 
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Figure 33. Performance of the Bundle with Simulation Heaters for 
Increasing Heat Flux in R-114 with 2% Oil 
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Figure 34. Performance of the Bundle with Simulation Heaters for 
Decreasing Heat Flux in R— 114 with 2% Oil 



81 



100000 




O O O 

O O t- 

O t- 

T— 

(2IAI /M) xnid 1V3H 

Figure 35. Performance of Tube 1 for Increasing/Decreasing Heat Flux in 
R-114 with 3% Oil 
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Figure 36. Performance of the Bundle with 
Increasing Heat Flux in R— 114 with 3% Oil 
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Figure 37. Performance of Tubes 1 and 2 for Decreasing Heat Flux in R- 
114 with 3% Oil 
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Figure 38. Performance of Tubes 1, 2, and 3 for Decreasing Heat Flux in 
R-114 with 3% Oil 
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Figure 3 9. Performance of Tubes 1. 2, 3, and 4 for Decreasing Heat Flux 
in R-114 with 3% Oil 
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Figure 40. Performance of All Five Tubes for Decreasing Heat Flux in R- 
114 with 3% Oil 
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Figure 41. Performance of All Five Tubes with Active Pairs for 
Decreasing Heat Flux in R— 114 with 3% Oil 
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Figure 42. Performance of the Bundle with Simulation Heaters for 
Decreasing Heat Flux in R-114 with 3% Oil 
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Figure 43. Comparison of Tests One to Seven for Tube 1 for Decreasing 
Heat Flux in R-114 with 3% Oil 
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Figure 44. Performance of Tube 1 for Increasing/Decreasing Heat Flux in 
R-114 with 6% Oil 
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Figure 45. Performance of the Bundle with Simulation Heaters 
Increasing Heat Flux in R-114 with 6% Oil 
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Figure 46. Performance of Tubes 1 and 2 for Decreasing Heat Flux in R- 
114 with 6% Oil 
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Figure 47. Performance of Tubes 1, 
R-114 with 6% Oil 



2, and 3 for Decreasing Heat Flux in 
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Figure 48. Performance of Tubes 1, 
in R-114 with 6% Oil 


2 , 3 , and A for 


Decreasing Heat Flux 
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Figure 49. Performance of All Five Tubes for Decreasing Heat Flux in R- 
114 with 6% Oil 
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Figure 50. Performance of All Five Tubes with Active Pairs for 
Decreasing Heat Flux in R-114 with 6% Oil 
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Figure 51. Performance of the Bundle with Simulation Heaters for 
Decreasing Heat Flux in R-114 with 6% Oil 
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Figure 52. Comparison of Tests One to Seven for Tube 1 for Decreasing 
Heat Flux in R-114 with 6% Oil 
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Figure 53. Performance of Tube 1 for Increasing/Decreasing Heat Flux in 
R-114 with 10% Oil 
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Figure 54. Performance of the Bundle with Simulation Heaters for 
Increasing Heat Flux in R-114 with 10% Oil 
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Figure 55. Performance of Tubes 1 and 2 for Decreasing Heat Flux in R- 
114 with 10% Oil 
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Figure 56. Performance of Tubes 1, 2, and 3 for Decreasing Heat Flux in 
R— 114 with 10% Oil 
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Figure 57. Performance of Tubes 1, 2, 3, and 4 for Decreasing Heat Flux 
in R-114 with 10% Oil 
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Figure 58. Performance of All Five Tubes for Decreasing Heat Flux in R- 
114 with 10% Oil 
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Figure 59. Performance of All Five Tubes with Active Pairs for 
Decreasing Heat Flux in R-114 with 10% Oil 
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Figure 60. 
Decreasing 



Performance of the Bundle with Simulation 
Heat Flux in R-114 with 10% Oil 
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Figure 61. Comparison of Tests One to Seven for Tube 1 for Decreasing 
Heat Flux in R-114 with 10% Oil 
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Figure 62. Comparison of Test One for Increasing Heat Flux in R-114 /Oil 
Mixtures 
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Figure 63. Comparison of Test One for Decreasing Heat Flux in R-114 /Oil 
Mixtures 
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Figure 64. Comparison of Tests One to Seven Tube One for Increasing Heat 
Flux in R-114/0il Mixtures 
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Figure 65. Comparison of Tests One to Seven Tube One for Decreasing Heat 
Flux in R— 114/Oil Mixtures 
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Figure 66. Mean Bundle Heat— Transfer Coefficient for Increasing Heat 
Flux in R-114/0il Mixtures 
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Figure 68. Test One Comparison of Turbo-B, Smooth, Finned, and High Flux 
Tube Bundles for Decreasing Heat Flux in Pure R-114 
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Figure 69. Ratio of Mean Bundle Heat-Transfer Coefficient With Oil to 
Heat-Transfer Coefficient Without Oil for Different Oil Percentages at 
a Heat Flux of 15 kW/m 2 
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Figure 70. Ratio of Mean Bundle Heat-Transfer Coefficient With Oil to 
Heat-Transfer Coefficient Without Oil for Different Oil Percentages at 
a Heat Flux of 30 kW/m 2 
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Figure 71. Ratio of Mean Bundle Heat-Transfer Coefficient of Enhanced 
Tube to Heat-Transfer Coefficient of Smooth Tube for Different Oil 
Percentages at a Heat Flux of 15 kW/m 2 
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Figure 72. Ratio of Mean Bundle Heat-Transfer Coefficient of Enhanced 
Tube to Heat-Transfer Coefficient of Smooth Tube for Different Oil 
Percentages at a Heat Flux of 30 kW/m 2 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

Nucleate boiling data of R-114 at atmospheric pressure were obtained 
using a small bundle of Turbo-B copper tubes. The data were obtained for 
both increasing and decreasing heat flux and at different oil 
concentrations. Based upon the results pertaining to this particular 
bundle and apparatus, the following conclusions may be made: 

1. Natural Convection Region 

a. For a single upper tube, a second lower tube directly below 
when turned on does increase the heat transfer performance of the upper 
tube, however when additional lower tubes are heated no net increase in 
performance occurs . 

b. The presence of heated lower tubes in the bundle reduces the 
incipient boiling point of the upper tubes and the tubes tend to nucleate 
'in order' (ie. top tube first, bottom tube last). 

c. The effect of adding oil to the refrigerant (up to 10%) 
reduces the heat-transfer coefficient slightly (approximately 10-15%) due 
to changes in the fluid properties. 

2. Boiling Region 

a. For pure R-114, the presence of heated lower tubes on the top 
tube causes no enhancement at high heat fluxes (> 20 kW/m 2 ) , but at low 
heat fluxes (< 20kW/m 2 ) , there is a significant enhancement due to 
convective effects . 
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b. At very low heat fluxes (<2 kW/m 2 ) , the presence of oil has 
little effect on the heat transfer performance of the top tube in the 
bundle. At higher heat fluxes (> 2kW/m 2 ) , the performance is enhanced by 
10-15% at low concentrations, but is degraded up to 20% at 10% oil 
concentration at the highest heat fluxes. 

c. At typical operating heat fluxes (15-30 kW/m 2 ) , the bundle 
performance is reduced between 5—15% with oil. 

B. RECOMMENDATIONS FOR FUTURE WORK 

1. Conduct experiments with varying pool height, but keep the local 
pressure at each tube constant by simul taneously varying the vapor 
pressure above the pool. 

2. Additional experiments with R-113 and R-114 should be conducted to 
investigate explosive (R-114) and partial (R-113) incipience at the onset 
of nucleation varying the time at the incipience. 

3. Some instrumentation should be added such that the flowrates 
through the bundle can be determined. From these measurements, vapor 
quality can be determined. 

4. Metal guide plates should be manufactured and placed on each side 
between the simulation tube bundle and the tube bundle itself. This 
further channels the flow of refrigerant thru the bundle at high heat 
fluxes . 

5. Attention needs to be given to the question of refrigerant 
disposal. There are reclamation projects undertaken by most 
manufacturers; however, a method still needs to be found to remove the 
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refrigerant from the apparatus into a container suitable for such 
reclamation . 

6. A high speed camera should be used to study the nucleation process 
and circulation patterns in more detail in the bundle. Neutrally buoyant 
particles might be placed in the pool to facilitate study of circulation 
patterns within the bundle. 
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APPENDIX A: LIST OF DATA FILE 



Table 3. DATA FILE NAMES FOR TURBO-B 
TUBE BUNDLE EXPERIMENTS 



FILE NAME 


NUMBER OF 
DATA POINTS 


NUMBER OF 
INSTRUMENTED TUBE 


PERCENT 
OF OIL 


NUMBER OF 

ACTIVE 

PAIRS 


NUMBER OF 

SIMULATION 

HEATERS 


TBI0001H 


34 


1 (15) 


0 


0 


0 


TBD0001H 


38 


1 (15) 


0 


0 


0 


TB I 0001 1 


25 


1 (13) 


0 


0 


0 


TBD0001 1 


22 


1 (13) 


0 


0 


0 


TB I 0001 J 


34 


5 


0 


0 


0 


TBD0001J 


25 


1 


0 


0 


0 


TBI0002A 


38 


2 


0 


0 


0 


TBD0002 


25 


2 


0 


0 


0 


TB I 0003A 


3t 


3 


0 


0 


0 


TBD0003 


27 


3 


0 


0 


0 


TBI 0004 


34 


4 


0 


0 


0 


TB00004 


24 


4 


0 


0 


0 


TBI0005 


34 


5 


0 


0 


0 


TBD0005 


24 


5 


0 


0 


0 


TBI 0006 


32 


5 


0 


5 


0 


1 TBD0006 


26 


5 


0 


5 


0 
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FILE NAME 


NUMBER OF 
DATA POINTS 


NUMBER OF 
INSTRUMENTED TUBE 


PERCENT 
OF OIL 


NUMBER OF 

ACTIVE 

PAIRS 


NUMBER OF 

SIMULATION 

HEATERS 


TBI0007 


31 


5 


0 


5 


5 


TBD0007 


24 


5 


0 


5 


5 


TBI 0101 


33 


1 


1 


0 


0 


TBD0101 


23 


1 


1 


0 


0 


TBI 0107 


29 


5 


1 


5 


5 


TBD0107 


26 


5 


1 


5 


5 


TBI0201 


38 


1 


2 


0 


0 


TBD0201 


22 


1 


2 


0 


0 


TBI 0207 


25 


5 


2 


5 


5 


TBD0207 


20 


5 


2 


5 


5 


TBI0301 


34 


1 


3 


0 


0 


TBD0301 


25 


1 


3 


0 


0 


TBD0302 


24 


2 


3 


0 


0 


TBD0303 


22 


3 


3 


0 


0 


TBD0304 


24 


4 


3 


0 


0 


TBD0305 


24 


5 


3 


0 


0 


TBD0306 


21 


5 


3 


5 


0 


TBI 0307 


26 


5 


3 


5 


5 


TBD0307 


21 


5 


3 


5 


5 


TBI0601 


32 


1 


6 


0 


0 
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FILE NAME 


NUMBER OF 


NUMBER OF 


PERCENT 


NUMBER OF 


NUMBER OF 




DATA POINTS 


INSTRUMENTED TUBE 


OF OIL 


ACTIVE 


SIMULATION 










PAIRS 


HEATERS 


TBD0601 


22 


1 


6 


0 


0 


TBD0602 


22 


2 


6 


0 


0 


TBD0603 


24 


3 


6 


0 


0 


TBD0604 


22 


4 


6 


0 


0 


TBD0605 


24 


5 


6 


0 


0 


TBD0606 


21 


5 


6 


5 


0 


TBI0607 


27 


5 


6 


5 


5 


TBD0607 


21 


5 


6 


5 


5 


TBI 1001 


33 


1 


10 


0 


0 


TBD1001 


22 


1 


10 


0 


0 


TBD1002 


27 


2 


10 


0 


0 


TBD1003 


24 


3 


10 


0 


0 


TBD1004 


24 


4 


10 


0 


0 


TBD1005 


23 


5 


10 


0 


0 


TBD1006 


16 


5 


10 


5 


0 


TBI 1007 


24 


5 


10 


5 


5 


TBD1007 


17 


5 


10 


5 


5 


TBI 0001 A 


24 


1 


0 


0 


0 


TBI0001C 


27 


1(15) 


0 


0 


0 


TBI 00010 


33 


1(13) 


0 


0 


0 
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FILE NAME 


NUMBER OF 
DATA POINTS 


NUMBER OF 
INSTRUMENTED TUBE 


PERCENT 
OF OIL 


NUMBER OF 

ACTIVE 

PAIRS 


NUMBER OF 

SIMULATION 

HEATERS 


TBI0001E 


33 


1 


0 


0 


0 


TBIOOOIF 


38 


1 


0 


0 


0 


TBIOOOIG 


38 


1 


0 


0 


0 
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APPENDIX B: SAMPLE CALCULATIONS 



Data set number 1 Tube 1 of experiment TBD1005 (Turbo-B tube, 

decreasing heat flux, 10% oil concentration, test 5) was used for the 
sample calculations in order to validate the program used for data 
acquisition DRP4RH. The working fluid was R-114. 

1. Test tube dimensions 

D tc = 11.60 mm 
D 0 “ 14.15 mm 
Di - 12.70 mm 
L = 203 . 2 mm 
Lu = 25.4 mm 

2. Measured Parameters 



T1 - 


10. 


.62 


°c 


T2 = 


10. 


.94 


°c 


T3 - 


9. 


.96 


°c 


T4 - 


10. 


.90 


°c 


T5 = 


10. 


.04 


°c 


T6 - 


9. 


.07 


°c 


Tldl 


f 


>.27 


°c 


Tld2 




>.21 


°c 


Aas «= 3 , 


.513 


V 
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Vas = 3.189 V 



3. Calculations 

The heaters power is first calculated for 
q = Vas(V)xAas(V)x60(V/V)xl(A/V) 

Note: The multiplication factors of volts and amp sensors are 60 and 
1, respectfully. 

Therefore : 

q * (3 . 189) (3. 513) (60V/V) (1A/V) 

q * 672.19 Watts 

The tube inside wall temperature is obtained from the average of all 
six thermocouple readings. 

6 

V1 6 A. D 

72-1 

- 1/6(10.62 + 10.94 + 9.96 + 10.90 + 10.04 + 9.07) 

- 10.25 °C 

The tube outside temperature is calculated by knowing the inside wall 
temperature using Fourier's Conduction Law. Uniform radial conduction is 
assumed . 

*[in(JB2.>] 

rp rp CC 

vo " w " 2* (* cu ) (L) 

where the second term on the right hand side is the Fourier conduction 
term. If we define this term as 
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q[ln(-£s-)] 

U tc 

2n (k^) (L) 



and 



* b - T uo - Tsat c 



where k cu is the thermal conductivity of copper and is calculated as 
follows 



k cu - 434.0 - [0.112 (T wj )] 



k cu - 434.0 - [0.112(283.25)] 



k cu - 4 02.28 W/mK 
now 



T ■- 

W1 



672.15 [In ( ^ 'gg ) 3 
2 n (4 02.28) (.2032) 



T vo - (10.25 - . 2601) ° C 



T uo - 9 . 9 8 c C 



The liquid saturation temperature at the top of the tube bundle is 
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Tsat 



tldl + Tld2 
2 



Tsa t 



2.27 + 2.21 
2 



Tsat - 2.24 °C 

In order to calculate the local saturation temperature for each tube, 
correction factors are needed to account for hydrostatic pressure 
differences between the tube locations and the liquid free surface. This 
difference is calculated by: 

AP - p (g) (ht) 

For the top tube in the bundle which is 0.124 m below the thermocouple 
measuring pool temperature. 

AP - 1523.12 (9.81) (0.124) 

A P - 1852 .78 Pa 

For 1852.78 Pa pressure difference, corrected saturation temperature is 
obtained by adding 0.04 °C (from standard tables for R-114) to Tsat. 
Corrected Tsat is: 

Tsat c - (2.24 + 0.04)°C 
Tsa t c - 2.28° C 
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Therefore, the wall superheat can be obtained by the following: 



- Two - Tsat c 
S b - (9.98 - 2.28) ° C 
6^ - 7 .70° C 

Now that the wall superheat is known, we need to calculate the heat 
flux and the heat-transfer coefficient. To do this, we know that the tube 
is 12 inches long and is heated in a eight inch center portion of the 
tube. The unheated lengths of the tube are a one inch and a three inch 
section on opposite ends of the tube. These unheated lengths have a fin 
effect during the heat transfer process to the evaporating refrigerant. 
In order to account for this, the following procedure was adopted for both 
one and three inch sections. Calculations are shown below for the one 
inch section. Heat transfer from the unheated end is calculated as heat 
from the base of the fin: 

Qt - Uh b ) (p) ( k CT ) (A c )]°-*(0 b ) ( tanh ( (n) (L c ) ] ) 

where 

P * 7T (D 0 ) 

- 7T (.01415) m 

- .04445 m 

now 

^ = n/U(D z 0 - D 2 j ) 

- 7t/4 ( . 01415 2 - . 01270 2 ) 

- 3.0578 x 10" 5 m 2 
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The corrected length of unenhanced surface at the end was calculated as 
follows 

Lc - Lu + (t/2) 

= 0.0254 + [(0. 01415-0. 0127)/2] 

- 0.0258 m 

h b is the natural convection heat transfer coefficient of the fin like 
ends and was calculated by using Churchill-Chu [Ref. 22] correlation for 
natural convection for a smooth cylinder, as modified by Pulido [Ref. 27], 



[ 



~ [ 0 . 6 + . 3 87 • 



[<7(P) (4>o) (9 b ) (tanh (nL c ) ) ] 
V (a) (L c ) ( n ) 



1/6 



[l+I- 



,559 

Pi 



9/16 



8/27 



where 



n 



(h b ) (p) 
(k cu ) (A c ) 



Therefore an iterative technique was necessary to calculate h^. The 
iterative technique used was to assume h^ was 190 W/m 2 K and continue the 
iteration until successive values are within 0.001 of each other. The 
fluid physical properties are calculated at the vapor mean film 
temperature, given by the following equation. 



T film 



Tsat c + TVo 
2 
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2.28 + 9.98 

1 film " ^ 



T /jJm - 6.13° C - 279 . 13 c K 



For R-114, the physical properties are given in the program by 
Dynamic viscosity . Tfjim in °K 

H = exp [-4 . 4636 + (1011 .47/T film ) ]xlCT 3 
H = 430.927xl0 -6 kg/m s 
Specific heat . Tfji m in °K 

C p = [0.40188 + 1.65007xl0" 3 (T film ) + 1 . 51494xlO- 6 (T f ilm 2 ) - 
6 . 67853xlO -10 (T f j[ m 3 ) ]xl0 3 
C p = 966.31 J/kgK 
Density . T filni in °K 

p - 16.0184533 (36.32 + 61 . 146414i|» 1/3 + 

17 . 47 6838<|r 1/2 + 1 . 119828>|i 2 ) 



where 



* - 1 



[1.8 (Tfilm) ] 
753.95 



and 



<| f - .333 

p - 1512.09 -^2 
m 3 

Thermal conductivity of R-114, Tfji m in °C 
k - 0.071 -(0.000261)(T film ) 



136 



k - 6 . 936xl0 -2 W/mK 



Prandtl Number 

Pr - [<Cp)p]/k 
Pr - 6.003 

Thermal Expansion Coefficient 

P - - (l/p) (Ap/AD 

p 279 . 03 - 1512.395 iS 
/n 3 

p 279 . 2J - 1511.824 iS 
m 3 

P - -(1/1512.395) [(.571) / (0.2)] 

P - 1 . 89x1 0 -3 (1/JO 
Kinematic viscosity 



v - JL 
P 

_ 430 . 927 Xl 0~ 6 
1512 . 09 

v - 2 . 8 49x1 0" 7 m 2 /s 



Thermal Diffusivity 
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a - 



k 

(p) C p 

6 .936 xiCT 2 
(1512.09)966.31 

a - 4 . 7 47xJ 0~ s m 2 /s 

Knowing the above properties, the heat-transfer coefficient h b , can be 
obtained by iteration 



h b = 362.57 W/m 2 K 



Knowing this we can calculate n 



n 



, (h b )( P ) °-5 

( k ^ (A c ) 



(362.57) (44 . 45 jcJ O' 3 ) j °- 5 
(402.28) ( 3 . 057 8x1 O' 6 ) 



n - 36.15 



then we can obtain qf 

q t - [ ( h b ) ip) (k cu ) ( A c ) ) °- 5 (6 b ) (tanh [ (n) ( L c ) ) ) 



q t - [(362.57) (.04445) (402.28) (3 . 057 8X10' 5 ) ] 05 
(7.70) (tanh [(36.196) (0.0258)]) 

q t - 2.51 W 



The corresponding results for the three inch section are 
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h b = 289.47 W/m 2 K 
q f - 1.24 W 



Therefore , the heat transfer through the heated length of the tube is 



q s — q — qf (1 inch section) — qf (3 inch section) 
q s = (672.19 - 1.24 - 2.51) W 
q s = 668.35 W 

and the heat flux and the heat transfer coefficient are as follows 

q" “ q s / A s 

- q s /((»)(Do)(L)) 

= (668. 35)/((ir)(0.01415)(. 2032)) 

= 7.398x10'' W/m 2 

and finally the heat transfer coefficient 



A t (T m " Tsat) 

h _ 668.35 

9 . 033x1 O' 3 (7 .70) 

h - 9 .605x1 0 3 W/m*K 
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APPENDIX C: UNCERTAINTY ANALYSIS 



The same data run (TBD1005) was chosen for the uncertainty analysis. 
Therefore, the measured and calculated parameters found in the sample 
calculation were used in this section. The unceratainty analysis 
performed was for a high heat flux, but the procedure could be performed 
at any heat flux to determine the uncertainty bands. All uncertainties 
are presented as a percentage of the calculated parameter. The 
uncertainty associated with the experimental parameters is calculated from 
the equation suggested by Kline and McClintock [Ref. 35]. For example: 

R « R(x 1> x 2f ,x n ) 

then 



6 R - I ( -#5- 8xj)M 

dx 1 1 



Sr 

dx y 



&x,) ; 



. + < 



dR 



t>*n> 



2T 0.5 



where 

SR = uncertainty of the desired dependant variable 
x n = measured variables 
6x n = uncertainty in measured variables 
The boiling heat-transfer coefficient is given by 



h 



_<?s 

A e (T vo -Tsat) 



where 



g[ln(-^>] 

rp T - UtC 

vo vi 2 n(k cu ){L) 
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In the above equation, the second terra on the right hand side is 
usually called the Fourier heat-transfer conduction term. If we define 



this as 



<y[ln(-^2-)] 

u tc 



2n(k cv ) (L) 



and 




With this notation, the uncertainty in the heat-transfer coefficient 
is obtained using the following equation. 



q= Vxl 

q- V(V)xI (V)x60 (V)xl (A/V) 
and the uncertainty is 

is . [<i i V.<ii,V ! 

q V I 

The accuracy in the voltage and current sensors are as follows 



bh 

h 



[(l2) 2 + (i^£)2 + (^k)2 + (M)2 + ( 

Q A s 0 b 0 b 



6Tsat ) 2 ] o.s 



where 
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6 Vas - ±0.05 V 



6Aas - ±0 . 025 A 



From the sample calculation section 

Vas = 3.189 V Aas = 3.513 V 

Therefore , 

_ [ ( AYsiJi) 2 + ( _&Aa£) 2 ] °-5 

q Vas ; Aas ; 

&Q _ r/ 0.05 \ 2 + / 0«025 x 2 iq,5 
q 3.189 3.513 

-^2 - 1.7 2 percent 
Q 



Calculation of the uncertainty of the surface area is as follows 
A s = tt(D 0 ) (L) 




L) rj L 



0.5 



Knowing the dimensions of the tube from the manufacturer and estimated 
inaccuraccies from work shop tools and human error, the unceratinty was 
calculated . 

Dimensions 

D 0 = 14.15 mm L = 203.2 mm 

Inaccuracies in measurements 
S D 0 = 0.1mm <5L = 0.2 mm 

Uncertainty analysis performed 
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[(4? £ ) 2+ <-^ : ) 2 ] oe 



. [(.. 0-1 ) 2 + (. 0-2 ) 2 ]0 . 5 
A* 14.15 203.2 J 



- - 0.7135 percent 

A s 

The uncertainty calculation for the Fourier conduction term given below 



... [ ( 6 ) 2 + ( ^ ^cu ^ 2 ^ ) 2 j 0 . 5 

4 > Q k cv L 

k cu was calcluated using 

k cu - 4 34.0- [ 0 . 112 ( T vi ) ] 

- 434.0- [0.112(283.25)] 
k cu - 4 02.28 W/mK 

and its uncertainty 



bk cu - [(0.112(ftr vi ) ) 2 ] 0 * 5 

<5T w j and <5Tsat are obtained using uncertainties in the thermocouple 
readings. Average wall inside temperature T W1 * was obtained taking the 
average of six thermocouple readings inside the tube wall. The 
uncertainty associated with this variable is 

Y6T . 2 

6 T wi - [€( vl ) ]°’ 5 
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where 6T gl * for each thermocouple was obtained by taking the difference 
between the measured wall temperature and the average wall temperature. 
Using this method, it has been attempted to try and take into acount some 
the uncertainty introduced by the fabrication procedure for the tube (ie. 
air gap) . For this particular heat flux the following 5T W1 * were found for 
tube thermocouples 1, 2, 3, 4, 5, and 6. 



5r 






0 

6 6 6 6 6 



0.5 



bT wi - .274° C 



The uncertainty level for all remaining thermocouple readings (ie. 
excluding those in the tube wall which was considered to have a higher 
uncertainty) was estimated to be +/— 0.1 °C corresponding to an emf of 
approximately 4 /iV. Saturation temperature was obtained by taking the 
average of two thermocouple readings and the uncertainty in this 
temperature was calculated from the following equation. 



* rp 2 

bTsat - [2(—^) ] 0 5 
bTsat - [2 (^l) 2 ] 0 - 5 
bTsat - 0.07 ° C 

Knowing the uncertainty in the temperatures, we can now calculate the 
uncertaines in the following: 
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6 k cu - [(0.112(6 T wi ) ) 2 ] 0 • 5 
bk cu - [(0. 112(283. 25) ) 2 ] °* 5 
bk cu - 31.724 W/mK 

Now we can calculate the uncertainty in the Fourier conduction term 



H ( & Q ) 2 , ( ^ ^ CV ^2, ( 8 ) 2 ] o . s 

4> q k cv L 






- 8.07 2 percent 

4 > 



from the sample calculations we know that 



4> - 0.267 9 ° C 6<J> - 0.2679 (0.08072) - 0.0216 ° C 



therefore : 



fr-k _ [( 6rsat ^ 2 j 0 .5 



A 



- [ (0 . 017 2 ) 2 + (0.07 135) 2 + {SLS±) 2 + ( °’°^ 6 ) ^ ) *1 0,5 

h 7.70 7.70 7.70 






7.42 percent 



Finally the calculation of wall superheat temperature 



- T vo - Tsat 



4 r ^ " [( Ar 2 ) 2 + (■ 6 ^ a£ ) 2 ] 0 - 5 



145 



&e * _ r , 0.274 / .07 , 2-10.5 

6 b 7.7 7.7 

86 , 

- - 3.67 percent 
0 *> 

Table 4 shows the results of the uncertainty analysis performed. The high 
and low heat flux correspond to the approximate values of 7.5x10* W/m 2 and 
lxlO 3 W/m 2 respectively. Note that the highest uncertainty (over 50%) is 
in the wall superheat at very low heat flux (1 kW/m 2 ) . This is due to the 
very low measured value of wall superheat (0.54 °C) which can not be 
accurately measured. Thus higher uncertainty occurs at very low heat 
fluxes. However, once the wall superheat gets higher (at higher heat 
fluxes) the uncertainty in wall superheat decreases significantly (to 
about 4%) indicative of the fact that the measure wall temperature is 
relatively more accurate. 
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Table 4. UNCERTAINTY ANALYSIS RESULTS 



VARIABLE 


HIGH HEAT FLUX 


LOW HEAT FLUX 


e* 


7.7 


0.54 


T wi 


10.25 


2.80 


^sat 


2.24 


2.20 


6V as/Vas 


1.57% 


14.6% 


6Aas/Aas 


0.712% 


7.9% 


Sq/q 


1.72% 


16.6% 


& D 0 /D 0 


0.707% 


0.707% 


S L/L 


0.098% 


0.098% 


5 As/ A^ 


0.714% 


0.714% 


Sk cu /k cu 


7.89% 


7.89% 


50^/e* 


3.67% 


52.37% 


6h/h 


7.42% 


28.93% 
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APPENDIX D: OPERATING PROCEDURE 



A . SYSTEM STARTUP 

1. Power to the 28 kW (8 ton) refrigeration unit is provided by the 
breakers located in the main distribution panel located in tzhe laboratory. 
These breakers were never secured. However, if power to this panel was 
lost, then these breakers must be reset. 

2. Turn the switch on the refrigeration unit control panel, located 
in front of the refrigeration unit to the "auto” position after passing 
through "on" position. This switch is always left on, unless unit was 
taken down for long repairs. 

3. Push the start button in the control box for the recirculation 
pump. This control box is located on the bulkhead above the recirculation 
pump in the outside area adjacent to the refrigeration unit. 

4. Set the desired temperature on the roughly graduated Fahrenheit 
scale on the control panel thermostat. It requires approximately one hour 
to chill the sump to - 15 °C. The thermometer in the ethylene glycol/water 
mixture (sump) must be monitored to ensure the desired sump temperature is 
attained and maintained. Slight adjustments in the refrigeration unit 
thermostat can be expected due to the coarseness of its scale. 

5. Energize the desired pumps by switching on the breakers in the main 
distribution power panel. Once the power is energized in the main power 
panel and after ensuring the pump suction valves are open, turn on the 
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pump motors by pushing down on the arm of the appropriate breaker box for 
the pumps located on the bulkhead next to the ethylene glycol/water sump. 
The pumps are marked "auxiliary condenser" (Pump #2) and "instrumented 
tube condenser" (Pump #1) t respectively in the breaker box. Flow in the 
auxiliary condensate coils can be controlled with the individual globe 
valves located at the coil penetrations on the apparatus. The auxiliary 
condenser coils will produce the fastest adjustments to the system 
pressure . 

6. Energize the heater variac(s) desired by switching on the breakers 
(Bundle, and Simulation (for test 7 only)) in the main distribution power 
panel and individual breakers for each of these in the power distribution 
box (near apparatus). 

7. After ensuring that the breakers for the heated tubes desired are 
in the "on" position, follow the experimental procedures for normal 
operation outlined in Chapter IV. 

B. SYSTEM SHUTDOWN 

1. Turn all variacs to the zero position and switch off all breakers 
in the power panels. 

2. If apparatus will not be operated for an extended period, turn the 
switch on the refrigeration control panel to the "off" position after 
passing through "on". 

3. Allow the recirculation pump to operate for at least five minutes 
after switching off the refrigeration pump unit to dissipate any back 
pressure in the system. 
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4. Turn the breakers for the pumps to the off position at the switch 



boxes, and then secure the power at the main distribution power panel. 



C . EMERGENCY SHUTDOWN 

1. Secure all power at the main distribution power panel 

2. Evacuate building 

3. Call Fire Department 
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APPENDIX E: PROGRAM DRP4RH 



1000 > 

10041 

1006' 

10:: ' 

1016 

10Z0 

1004 

10Z3I 

103: 
1026 
1040 
1044 
1048 
1052 
1056 
1060 
1064 
10681 
10721 
1076 
1080 
1084 
1088 
1092 
1096 
1100 
1104 
1108 
1112 
11161 
1120 
11241 
1128 
1132 
Ldt C ( 
1136 
11401 
1144 I 
1148 
1152 
1156 
1 160 * 
1164! 
1168 
1172 
1176 
1 1 80 * 
1184 
1188 
1192 
1196 
1200 
1204 
1208 
1212 
1216* 
1220 ' 
1224 
1228 



FILE NAPE DRP4RH 

DATE November 22, 1388 

REVISED FEB 1992 (R. HAAS) 



BEEP 
PRINTER 
I dp-0 



IS 1 



PRINT USING 
PRINT USING 



‘ 4 X , 
‘ 6x , 



PRINT USING " 6 X , " 
PRINT USING M 6 X , “ 
PRINT USING " 6 X , " 
PRINT USIN6 "6X, " 
PRINT USING " 6X , " 
PRINT USING " 6X , " 
PRINT USING " 6X , * 



"Select option default is 0 • '* “ ** 

“0 Taking data or re-proce ss 1 ng previous data" 

" 1 Plotting data on Leg-Log 

''2 Plotting data on Linear**"" 

M 3 Purge 
“4 Fixup""" 

"5 Move""" 

"6 Comb""" 

"7 Read Plot""" 



IDP IS A PROGRAM VARIA8LE TO SELECT A SU8R0UTINE 
INPUT I dp 

IF Idp-0 THEN CALL Main 
IF I dp- 1 THEN CALL Plot 
IF Idp-2 THEN CALL Plin 
IF Idp«3 THEN CALL Purg 
IF Idp-4 THEN CALL Fixup 
IF Idp-5 THEN CALL Move 
IF Idp-6 THEN CALL Comp 
IF I dp*7 THEN CALL Readplot 
END 

SU8 Main 

ICAL-THERMOCOUPLE CALI8RATI0N 
COM /Cc/ C( 7 ) 

DIM Em f < 35 > , T( 35 ) , D 1 a( 6 ),D2a<6 ),Dia(6 ) , Doa( 6 ) , La( 6 ),Lua(6 >,Kcua< 6 ) ,Et ( 19 ), 
20 ) , Uo 1 1 ( 2 > , Amp ( 1 1 ),Tua(5),Tu(5),Theta(5),Thetab(5),Q(5),Ql<5),Qdp(5 ) 

DIM Ht ube ( 5 > , Tn< 5 > , Tp < 6 ) 

THERMOCOUPLE ARRAY (C( )) INITIALIZATION 

DATA 0.10086091,25727.94369,-767345.9295,78025595.81 
DATA -9247486589, 6. 97688E+ 11 ,-2.661 92E+ 13,3. 94078E + 1 4 
READ C< ♦ ) 

PRINT HEAOER AND INITIALIZE TIME CLOCK 

PRINTER IS 701 

BEEP 

INPUT "ENTER MONTH, DATE AND TIME < MM DD ; HH • MM SS ) " , Da t e$ 

OUTPUT DIRECTEO TO DATA AQUISITION SYSTEM (HP 3497A) 

OUTPUT 709, "TD" .Dates 
OUTPUT 709. "TD" 

ENTER 709 : Dat e$ 

PRINT 

PRINT " Month, date and time -";DateS 

PRINT 

PRINT USING ** 1 0X , ** "NOTE Program name DRP4RH" " " 

8EEP 

DN IS THE VAPIA8LE FOR DISC NUM8EF FOR RECORD KEEPING ONLY 

INPUT “ENTER DISK NUM8ER " , Dn 

PRINT USING " 1 6 X , ** ** D 1 s k number - "".ZZ-.On 
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12Z-2 BEEF 
1236 In-0 

124C INPUT "ENTER INPUT MOOE <0- 34 S7 A, 1-FILE ) 0-OEFAULT ",1m 
1244 ' 

1248' INPUT MOOE ZERO IS FROM THE OATH AOUISITION SYSTEM 
12S2 IF In-0 THEN 
1256 8EEP 

1260 INPUT "GIVE A NAME FOR THE RAW OATA F I LE “ , 02 f i 1 eS 

1264 PRINT USING " 1 6 X , “"File name " \14A" . 02f : leS 

1268' 

1272' CREATE 80AT FILE ON THE MASS STORAGE MEOIA 

1276 CREATE GOAT O2filef,60 

1280' CREATE AN INPUT/OUTPUT LINK TO OPEN FILES 

1284 ASSIGN e p ile2 TO D2file$ 

1288' 

12921 CREATE OUMMY FILE UNTIL Nr u n KNOWN 

1296 D1 f l le$- "OUMMY" 

1300 CREATE BOAT DlfileS.60 

1304 ASSIGN ©Filel TO Dlfile* 

1308 OUTPUT ©Filel ;0eteS 

1312' 

1316! CREATE A PLOT FILE 

1320 BEEP 

1324 INPUT "GIVE A NAME FOR THE PLOT FILE",PfileS 

1328 CREATE 80AT PfileS,30 

1332 ASSIGN ©Plot TO PfileS 

1336 BEEP 

1340' 

1344 ' IOTC - NUMBER (TOTAL) OF DEFECTIVE THERMOCOUPLES 

1348 INPUT "ENTER NUM8ER OF OEFECTIVE TCS ( 0-DEFAULT >', I dtc 

1352' LDTC - LOCATION OF DEFECTIVE THERMOCOUPLE 

1356' 

1360 IF Idtc-0 THEN 

1364 PRINT USING "l6X,"No defective TCs exist"" 

1368 ELSE 

1372 PRINT USING " 16X , " "Oef ect l ve Thermocouples Indicated by -99.99"* 

1376 ENO IF 

13801 

1384 BEE P 

1388' DEFECTIVE THERMOCOUPLES MAY BE IN CHANNELS 40-69 

1392' THEMOCOUPLES ARE ENTEREO AS OEFECTIVE BY COMPUTER CHANNEL NO. 

1396! JDT C-COUNTER IN LOOP FOR OEFECTIVE THERMOCOUPLES 

1400' 

1404 IF I dtc >0 THEN 

1408 FOR Jdtc-0 TO Idtc-1 

1412 INPUT "ENTER OEFECTIVE TC LOCATION (8Y COMPUTER CHANNEL NUMBER ) 

" , Ldtc( Jdt c ) 

1416 8EEP 

1420 NEXT Jdtc 

1424 ENO IF 

1428 PRINTER IS 701 

1432 OUTPUT ©Filel ;Ldtc(* ) 

1436 ' 

1440 ! In- 1 option (THIS OPTION ALLOWS OATA ENTRY WITH OATA PILE' 

1444 ELSE 
1448 BEEP 

1452 INPUT "GIVE THE NAME OF THE EXISTING OATA F I LE " , D2 f : i eS 

1456 PRINT USING * 1 SX , " "F i le name " " , 14A" ; 02 f l let 

1460 ASSIGN ©file2 TO D2fiie* 

1464 ENTER @File2.Nrun 
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1468 ENTER ®F i 1 e2 . 2ol , LJtc 1 •), I t t , Boo , Nht , Na t d , N r t , Zorr 

14“’: SEEP 

14^5 INPUT “GIVE A NAME FOR PLOT FILE” ,Pf i leS 

1490 CREATE BOAT Pfile$,30 

U64 ASSIGN SPlot TO P'lie* 

1408 PRINT USING * 1 6 < , “ “ Th i 3 data set taken on ” " , l 4 A ” . Do l dS 

1482 SEEP 

I486 PRINTER IS l 

1500 PRINT USING -4* "SELECT TU6E TYPE 

1504 PRINT USING "6X,""0 SMOOTH" “ - 

1508 PRINT USING ” 6 X , " “ 1 FlNNEDt 19/IN ) 

1512 PRINT USING “6K,"”2 HIGH FLUX 

ISIS PRINT USING ”6X,""3 TURBO-B 

1520 INPUT Itt 

1524 END IF 

1528 IF In-1 THEN GOTO 1768 
1532 PRINTER IS 1 
1536 ' 

1540 IF In-0 THEN 



1544 


PRINT USING 


" 4X , " " 


Select tube type""" 


154 6 


PRINT USING 


*6X, " " 


0 


Snooth """ 


1552 


PRINT USING 


" 6X , "" 


1 


FINNED 19/IN (DEFAULT ) 


1556 


PRINT USING 


"6X , ” " 


■> 


HI6H FLUX""" 


1560 


PRINT USING 


"6X , " " 


3 


TURBO-B" “ " 


1564 


PRINT USING 


"6X,““ 


4 


GROWTH” " " 


1553 


PRINT USING 


"6X , " " 


5 


GROWTH" " " 


1572 


PRINT USING 


“ 6X , " " 


6 


GROWTH”" " 


1576! 


ITT-TUBE TYPE 






1580 


INPUT Itt 








1 5B4 


OUTPUT ©Filel ; Itt 







1 58B END IF 

1592 PRINTER IS 701 

1596 f Itt-2 

1600 PRINT USING " 16X , " “ Tube Type "",0D";Itt 

1604 1 

1 60B BEEP 
1612 Bop-0 

1616 INPUT "ENTER BULK OIL X ( DEFAULT -0X ) ",Bop 
1620 OUTPUT ®F 1 1 e 1 ; Bop 

1624 PRINT USING "16X,""Bulk 01 1 X* " " , DD " . Bop 
1628! 

1632 BEEP 

1633 Ipo-1 

1634 INPUT "ENTER POOL HEIGHT ABOUE TOP TUBE (0-LESS THEN 5cn, i-5cn OR 6REATER 
(DEFAULT)", Ipo 

1635 OUTPUT ®F 1 1 e 1 , Ipo 

1636 PRINT USING "16X,""Pool he 1 ght - " " , DD " . I po 

1637 BEEP 

1639' NHT -NUMBER OF HEATED TUBES 
1640 Nht-5 

1644 INPUT "Enter nunper of heated 1 na trunen ted tubea( defaul t-5 )" ,Nht 
1 64B OUTPUT @Fi lei .Nht 

1652 PRINT USING ” I6X Number of heated instrumented tube s- " " , DD “ . Nh t 

1653 BEEP 

1655 Ipos-1 

1656 IF Nht- 1 THEN INPUT "WHICH POSITION IS THIS SINGLE TU6E ( 1 TO 5 DEF=l)",I 
poa 

1657 Ipos»Ipo5-l 

1658 BEEP 
1660' 
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166*1 Natp-Numbe^ of active ounny pai r 5 
1666 Natp-0 

1672 I N D UT "Ente^ nuwbe r of active Ounny pairs i De f au 1 1 -0 ) " , Na t p 
1E7E OUTPUT ®F 1 1 e 1 ; Na tp 

168C PRINT USING " 1 6X , " "Number of active dummy pa i rs - " “ , DC “ . Na tp 

16 84 BEEP 

1668' 

1692 * NRT-NUM6EF OP AOOED HEATEO TUBES TO ENHANCE bundle EFFECT 
1E9E N-t*0 

1700 INPUT “Enter number of added heated tubes from simulation heaters( Default- 
er , Nr t 

1704 OUTPUT @F i le 1 ; Nrt 

1708 PRINT USINE * 1 6X , “ " Number of added heated tudestfrom simulation heaters)*’ 
" ,00" ;Nrt 
1712 BEEP 
1716 ' 

1720' CORR IS CORRECTION FDR INSTRUMENTED TU8E HEI6HT 
1724 Corr-0 

1728 INPUT “WANT TO CORRECT TSAT FDR TUBE HEI6HT ( 0-YES( DEFAULT ), 1 -NO )*, Corn 
1732 IF Corr-0 THEN PRINT USINE “1EX,”TSAT is corrected instrumented heat 

ed tube heiQht*** 

173E IF Corr*l THEN PRINT USIN6 *1EX,**TSAT is NOT corrected for lnstrumen 

ted heated tube heiQht*** 

1740 OUTPUT GFileliCorr 
1744 BEEP 

1748' ILOU-INPUT MODE LIQUID, UAP0R,DR LIQUID UAPDR AVERAGE 
1752 Ilqv-0 

175E INPUT * SELECT ( 0-LIG( default ) , 1-UAP,2*( LIQ+UAP )/2 )* , I Iqv 
17601 

1764' DlA-Diameter at thermocouple positions (meters) 

1 7E8 DATA . 01 22 , 0 . 0098 , 0 . 0 10E , 0 . 0 1 1 E , 0 , 0 , 0 
1772 REAO Diet • ) 

177E Dl-0la( Itt ) 

1780' 

1784' D2*Diameter to base of fins (outside dia for smoot h )( meters > 

1788 DATA .0158,0.0125,0.0158,0.01415,0,0,0 
1792 READ 02a( * ) 

1796 D2-D2a( Itt) 

1800 ' 

1804' Di*In*ide diameter of unenhanced ends (meters) 

1806 OAT A .0132,0.0109,0.0116,0.0127,0,0,0 
1812 READ Dia( * ) 

1816 Di-Dia(Itt) 

18201 

1824' Do-Outside diameter of unenhanced ends (meters) 

1828 OATA . 01 5875, 0.01 25, 0.0 15875. 0.01 4 15, C,e,0 
1632 REAO 0oe(* ) 

1636 Oo-Doa(Itt) 

1840 ' 

1844' L-Length of enhanced surface (meters) 

184 8 DATA .2032, .2032, .2032, .2032, .2032, .2032, .2022 
1852 REAO Let • ) 

1856 L-La(Itt) 

I860' 

1664' Lu- CORRECTED Length of unsnhanceb su r face at the ends (METERS 

1868' LU-LFIN - THICKNESS/2 

1672 OATA .0261, .0254, .C2E4. 0.0258,0,0,0 

1878 READ Lua< * ) 

168G Lu=Lua(Itt) 

1864 i 
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1 8 8 B 1 LV-corrected length o f 2 inch finned lile end 
163: 01* LVfl! 6 > 

1696 DATA . 0765 , . 0762 , . CT72 , 0 . 0766 , 0 . 0 ,0 
1900 READ Lvfli * ) 

1904 Lv-Lva(Itt) 

1906 1 Kcue-Thermel Conductivit> of tube 
1912* data 401 ,G,e,e,c,c,0 
1916 ' READ h cue * ' 

1920' hcu-Kcua< Itt ) 

1924 A-P I ♦ ( Dc -2-Di'2 )/4 

1926 P-PI*Do 

1932 J- 1 

1936 Sx-0 

1940 Sy-0 

1944 S x s * 0 

1948 S x y *0 

1952 Repeat: I 

19561 

1960 IF In-0 THEN 

1964' Dtld-desired temperature of liquid 

19681 Dtld-47.5 iR-113 

1969 Dtld-2.2 *R-114 

1972 Ido-2 

1976 ON KEY 0,15 RECOVER 1952 

1980 PRINTER IS 1 

1984 PRINT USING " 4X SELECT OPTION " " " 

1988 PRINT USING '6X,“0«TAK6 DATA* “ " 

1992 PRINT USING "6X , “ " 1 -SET HEAT FLUX"’ 

1996 PRINT USING "6X , " “2-SET Tsat (DEFAULT SET FOP R-114)"" 

2000 PRINT USING “4X , “ “NOTE KEY 0 - ESCAPE"' - 

20041 Ido-desired option 

2008 BEEP 

2012 INPUT Ido 

20161 

2020 BEEP 

20241 s e t default value for input 

2026 IF Idc>2 THEN Ido-2 

2032 1 Take data option 

2035 IF I do-0 THEN 2440 

2040! 

20441 LOOP TO SET HEAT FLUX (FOR TOP INSTRUMENTED TUBE) 



2048 

2052 

2056 


IF Ido-1 THEN 
Dqdp-100000 

PRINT USING 


" 4X , " 


‘Odp 


QDPs l m 


Nr t 


Gdpaux 


C t o t ‘ " 
2060 


PRINT USING 


“ 4X , " 


“< W/m"2 ) 


( W/m'2 ) 




( LJ / r ~ 2 


( LI ) " “ “ 
2064 


Err* 1 













2066* Reset, read channel 25-30 , au t oma t i c scaling 

2072 1 Channel 25 = au> emps,26-sir amps,27«inst vcl t s , 26-s ir. volts, 29-au/ 

vol t s , 30-34= i ns t amps 

2076 OUTPUT 709; "AR AF25 AL34 UR5 " 

2080 FOR 1-10 TO 11 

2064 OUTPUT 705, 'AS SA' 

2066 ENTEF 709 . Amp t I ) 

2092 NEXT I 

2096 FOF 1=0 TC 2 

2100 OUTPUT 705- "AS SA M 

2104 ENTEF 709 . Uo 1 t ( I > 

2106 NEXT I 
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:ics 

2116 

2120 ' 

2124 

212E 

2152 

213E 

2140 

2144 

214E 

2152 

2156 

2160 

2164 

2168 

2172 

2176 

2180 

2184 

2188 

2192 

( 0 ) , Qdps i 
2196 
2200 
2204 
2208 
22 1 2 1 
2216 ' 

2220 

2224 

2226 

2232' 

2233 

2236 

2240 

2244 

2246 

2252 

2256 

2260 

2264 

T va t 

226e 

2272 1 

2276 

2230' 

rnocoup 1 € 

2264 

2 2 BE 

2292 

2296 

2300 

2304 

2306 

2312 

2 3 1 E 

2320 

2324 i 



FOR I*C TC 4 
OUTPUT 709. “AE SA*’ 

ENTEF 70S. Amp< 1 ) 

NEXT I 

Calculate actual heat flu> 

0(0 >*6C*Uolt(0 >*Amp( I p o s ) 

Oap(C)-O(0)/(PI*O2*L) 

Gsim*60*20*Uolt(l >*Amp ( 1 1 ) 

Qdpsim-Qsim/(FI*D2*.2032*5^ 

Qau*“60*20*Uolt(2 >«Amp( 10 ) 

Gdpau»*Gaux/(Pl* .0160* . 1 778 • 4 ) 

Otot*Q( 0 )«Nht+Gsim+Gau* 

Nrt-Gdps im/Qdp< 0 ) 

IF ABS( Aqdp-Oqdp ) >Err THEN 
IF Aqdp>Dqdp THEN 
BEEP 4000, .2 

ELSE 

8EEP 250, .2 
END IF 

IF Nr t < . 1 THEN Nrt-0 
IF Qdpaux<100 THEN Qdpaux-0 
IF Qdpsim<100 THEN Qdpsim*0 

PRINT USING "4X ,2<MZ .30E,2X >,2X , (MDO.DO ) , 2X ,2<ttZ .3DE,2X )“ ;Qdp 
m , Nr t , Qdpau* , Otot 
WA IT 2 
60T0 2076 
ENO IF 
END IF 

LOOP TO SET Tsat 
IF I do*2 THEN 

IF Ikdt-1 THEN 2240 
8 EEP 

INPUT ‘ENTER DESIRED Tsat < DEFAULT-47 . 5 C - R-113 >" ,0t Id 

INPUT “ENTER DESIRED Tsat (DEFAULT*:. 2 C - R- 1 14 ) " , Dt Id 

Ikdt-1 

01 dl *0 

01 d2-0 

Nn* 1 

Nrs-Nn MOD 15 
Nn*Nn+ 1 
IF N-s-1 THEN 

PRINT USING "4X,““ DTsat Tldl Tld2 Tlbb 

Tvab T la v 

END IF 

Read thermocouple voltages for vapor, liquid 
OUTPUT 705 . * AR AF 0 AL5 VR5 " 

Sanple each thermocouple 20 times and report temp for each the 
, vapor-0,1,2; liquid*3&4 
FOR >0 TC E 

S u m » C 

OUTPUT 705;“ AS SA“ 

FOR Ji*i TO 20 

ENTER 709 . E 1 i q 

Sum= Sum + E 1 i q 
NEXT Ji 
Emf(I ? - S u m / 2 0 
T ( I ) ■ F N T vsv( Emf ( I ) > 

NEXT I 

Compute average temperature of liquid 
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;:34 
2 j jE 
2342 
2344 
2346 
2352 
2356 
2360 
2364 
2368 
2 j 72 
2376 
2360 
2364 
2388 
2392 
2396 
2400 
2404 
2408 
2412 
2416 



T lav-l Ti 2 >*T( 4 >)• .5 

Compute average, temperature of vapo'- 

Tvav ] « ( 7 ( C ) + T ( 1 ) )/2 

Tvav2-T( 2 ) 

Tvav«( T(0 >♦ T ( 1 )+ T ( 2 > )/3 
IF ABS( T 1 av-Dt 1 d . 2 THEN 
IF TlavlDtlc THEN 
8EEP 4000, .2 

ELSE 

8 EEP 250, .2 
END IF 

ELSE 

IF A8S ( T 1 a v-Ot 1 d ) .1 THEN 
IF At ld>0t Id THEN 
EEEP 3000 ,.2 

E L S E 

BEEP 800, .2 
ENO IF 
ENO IF 
END IF 

Err 1-T lev-Dl dl 
Oldl-Tlav 
Err 2“ T vav-Dl d2 
D 1 d2*T vav 

PRINT USING * 4X , 7( NDDD . 00 , 3X >";Dtid,T(3),T(4 ),T(5),Tvavl ,Tvav2,Tla 



V 

2420 


WAIT 2 


2424 


GOTO 2252 


2428 


END IF 


2432 1 




2436 ' 


TAKE DATA IF In-0 LDDP 


2440 


IF Ikol-1 THEN 2452 


2444 


8EEF 


2448 


Ikol-1 


2452 


DUTPUT 709 i ’ AR AF0 AL5 UR5“ 


2456 


FOR 1-0 TD 5 


2460 


DUTPUT 709. “AS SA“ 


2464 


Sun-0 


2468 


FOR Ji-1 TO 20 


2472 


ENTER 709. E 


2476 


Sun-Sun+E 


2480 


IF I >2 THEN Et< Ji-i )-E 


2484 


NEXT Ji 


2488 


Kdl-0 


2492 


IF I>2 THEN 


2496 


Eave-Sun/20 


2500 


Sun-0 . 


2504 


FOP Jf-0 TO 19 


2508 


IF A8S( Et ( JI )-Eave X5.0E-E THEN 


2512 


Sun-Sun+Et ( Jk ) 


2516 


ELSE 


2520 


Kdl-Kril+1 


2524 


END IF 


2528 


NEXT 


2532 


IF 1.2 THEN PRINT USING “AX/^Kdl - “ “ ,DD*‘.kd] 


2535 


IF \ dl 10 THEN 


2540 


EEEP 


2544 


EEEP 


254e 


PRINT USING “4X,““Tcc nuen scattering in 


peat aete 


s e t “ “ “ 



dote 



re 
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rss: 

2556 
2560 
2564 
2566 
2572 
2576 
2580 
2584 
2588 
2592 
2596 
2600 
2604 
2606 
2612' 
2616 * 
2620 
2624 
2628 
2632 
2636 
26401 
2644 i 
2648 
2652 
2656 
2660 
2664 
2668' 
2672 
2676 
268C 
2684 
2688 
2692 
2696 
2700 
2704 i 
2708 ' 
2712 
2716 
2720 
2724 
272e 
2732 
2736 
2740 
2744 i 

2748 

2749 

2750 

2751 

2752 

2754 

2755 

2756 

2757 
2755 
27Ee 



GOTO 1 5BC 
END IF 
END IF 

Ewf ' I » = Sum/ \ 20-hdl ) 

NEXT I 

OUTPUT 709 . " AR AF 40 AL69 UR5" 

FOP 1-6 TO 35 

OUTPUT 709. “AE SA" 

Sum-0 

FOF J i-l TO 5 
ENTEF 705. E 
Sum-Sum+E 
NEXT Ji 
Emf l I )-5un/5 
NEXT I 

READ U0LTA6ES ( 27- I nst , 26-Sim , 29-Au> ) 

OUTPUT 709 ■ * AR AF27 AL25 UR5" 

FOR 1-0 TO 2 

OUTPUT 709. "AS SA " 

ENTER 709 ; Uol t ( I ) 

NEXT I 

READ CURRENTS (30-34-Inst tubes ; 3S-39-ACT IUE Dummy) 
OUTPUT 70S . " AR AF30 AL39 UR5“ 

FOR 1-0 TO 9 

OUTPUT 709. "AS SA" 

ENTER 70S . Amp ( I ) 

NEXT I 

Read Current s ( 25”Aux amps,26-$im amps) 

OUTPUT 709; "AR AF25 AL26 UR5 " 

FOR 1-10 TO 11 

OUTPUT 709; "AS SA" 

ENTER 70S ; Amp ( I ) 

NEXT I 
ELSE 

ENTER Ie2 ;Emf< • ),Uolt< ♦ ),Amp( * ) 

END IF 

CONUERT EMF’S TO TEMP , UOLT , CURRENT 
FOR 1-0 TG 35 

T ( I )-FNTvsv(Emf< I ) ) 

IF I >5 AND I dt c>0 THEN 
FOR I i-0 TO Idtc-1 

IF Ldt c ( Ii )* I -5 + 39 THEN T(I> — 93.99 
NEXT Ii 
ENO IF 
NEXT I 

Ntc-nr of thermocouples 
Nt c-6 

IP I pos >0 THEN 

C( I pos >*50»Uol t ( 0 )*Amp( I pos ) 

Twa^ Ipos '-0 

Jj-0 

Nd t c -0 

FOR I-l TO Nt c 

No* Ipos*6 + 5 + J j 
J J = J j «• 1 

I c AES( T ( Nn ) ) 95 THEN 
T ( Nn >--SS . 95 
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2761 


Ndtc-Ndic+1 


2762 


ELSE 


2763 


Twe( Ipos )=Twa< Ipos > + T<Nn> 


2764 


END IF 


2765 


NEY- I 


2766 


Twe( Ipos >-Tua< Ipos )/( 6-Ndtc ) 


2767 


GOTO 2620 


276S 


END IF 


2765 


FOR Ii-C TC 4 


2770 


Q - Ii >*60*Volt(C>*Anp(Ii ) 


2771 ' 


Twa-AveraQe temperature of the wall 


2772 


Twa ( 1 1 >“0 


2773 


Ndtc-0 


2774 


FOR 1*1 TO Ntc 


2776' 


Nn is counter in temp a-ray , ste*-t 


ermocouple in the tube ban! ) 


2780 


Nn- I i «6+ I +5 


2784 


IF AB5( T( Nn ) )> 99 THEN 


2788 


T( Nn >*-99 . 99 


2792 


Ndtc*Ndtc+l 


2796 


ELSE 


2800 


Twa( li >-Twa( Ii ) + T( Nn ) 


2804 


ENO IF 


2808 


NEXT 1 


2812 


Twa( Ii )-Twa( Ii >/( 6-Ndtc > 


2816 


NEXT Ii 


2820 


T I a v* ( T ( 3 >+T ( 4 ) )/2 


2821 


T vav*T ( 2 ) 


2824' 

2828' 


T va v* ( T ( 0 )+T ( 1 )+T< 2 ) )/5 


2829 ' 


T I a v*T ( 5 ) 


2832 


T cu*Twa( 0 > 


2836 


Kcu*FNKcu< Tcu > ITHERMAL CONDUCTIVITY OF 


2840' 


I IF CURVE FIT NOT AVAIL 



at 6 (this is the first th 



28441 FOURIER CONDUCTION EQUATION WITH CONTACT RESISTANCE NEGLECTED 
2848 FOR 1*0 TO 4 

2852 Tut I )* Twa ( I )-Q< I >*LOG( 02/01 >/ < 2*P I *Kcu*L ) 

2856 IF llqv*0 THEN Texs*Tlav 

2860 IF I Iq v* 1 THEN Texs-Tvav 

2864 IF IIqv-2 THEN Tex s- ( T I av + T ( 2 ) >* . 5 

2868 IF Corr-I THEN The t ab ( I )*Tw< I >-Tex s 

2872 IF Corr-0 THEN Thetab' I >*Tw< I >-( Tex.s* . 056+1 * . 129 > 'R-114 

2876 ' IF Co rr*0 ANO Ipo*l THEN Thetabt I >*Tu< I >-< Texs + .054 + 1* . 144 ) 



iR-Il 



2877 i 
113 
2880 
2884* 
2888* 
2892 
2896 
2900 
2904 
2908 
2912 
2516 
2920 
2924 
2928* 
2932 1 



IF Corr*0 AND Ipo-0 THEN Thetab( I )*Tu( I WTexs-l.C78+. 147*1 ) 



IR- 



NEXT I 



COMPUTE VARIOUS PROPERTIES 



TfiIr'“(Tu(0)-*-Tex5 )• . 5 
Rho**RNRho( T f l I n ) 
Mu*FNMu< Tf l lm ) 
f-RNK(Tfilm) 

Cp*FNCp ( T f l I m ) 
Beta^RNBetatTfilm ) 
Ni*Mu/Phc 
AIphs-K/< Pho*Cp ) 

Pr*Ni /A I pha 



i FILM TEflPERATUPE 
* DENSITY 
■VISCOSITY 

I thermal conductivity 

'SPECIFIC HEAT 
'THERMAL EXPANSION 
i v INEMATIC VISCOSITY 
'THERMAL DIFFUSIVITY 
IPRANDTL 



COMPUTE NATURAL-CONVECTIVE HEAT-TRANSFER coefficient 
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2936 1 F Of- UNENHANCEC END( S ) 

29^0 Lu'lus v 1 t t > 

2544 Hbar=19e 

2949 Fe c (Hbar*P/(hcu*A))'.5*Lu 

2952 7 anh-FNT anh ( Fe > 

2956 Theta( Ipos )-Thetab( Ipos )*Tanh-'Fe 

2980 Xx-(9.6l*Beta*Thetab( Ipos MOc'D^Tanh/CFe^Ni^Alpha) ' ' .18688? 

2 964 Vy-( l + < .559/Pr >- ( 9 / 16 ) > * ( 8/2? ) 

2 96 6 Hb arc-k/Do*( ,6 + .367*Xx/Yy)' , 2 

2972 IF AeS< (HDa--Hba^c )/Hbarc K .001 THEN 

2976 Hbar-( Hbar + Hbarc )♦ . 5 

2980 60TC 2946 

2984 ENO IF 



2 988 5 

2992' COMPUTE HEAT LOSS RATE THR0U6H UNENHANCED END(S) 

2996 01 ( 0 >-( Thetab( Ipos >*Tanh )• ( ( Hbar *P*Kcu* A )' . 5 ) 

3000 Qq-Q 1 ( 0 ) + 0q 

3004 2*2+1 

3006 IF 2-1 THEN 

3012 Lu-Lv 

3016 GOTO 2944 

3020 END IF 

3024 2-0 

3028 Qlpct-Qq/Q( Ipos > 

3032 Qq*0 

3036 As-PI *D2*L 

3040 FDR I i-0 TD 4 

3044 01( Ii )*Qlpct • Q( 1 1 ) 

304 8 Odp(Ii )-(0( Ii >-Ql< Ii ) )/As 

3052 Htube(Ii)*Qdp(Ii)/Thetab(Ii) 

3056 NEXT Ii 

3060 PRINTER IS 701 

3064 1 

3068' RECORD TIME DF DATA TAKING 

3072 IF In-0 THEN 

3076 OUTPUT 709, "TO" 

3080 ENTER 709.Told$ 

3084 END IF 

3085' CHURCHILl/CHU CORRELATION FDR NATURAL CONNECTION REGION 



3087 

3086 

3089 

3090 

3091 I 

3092 i 



Raa-9.8l*5eta*Thetab( Ipos )*(02 K3*Rno/(Mu*Alpha ) 
Oenom=< 1 . + ( .559/Pr )“( 9/ 16 ) r < 16/9 ) 

Nuc h* ( . 6+ . 387* ( Raa /Denom )"(l/6))*2 
Qch-K«Nuch*Thetab( Ipos )/<D2 > 

OUTPUT DATA TO PRINTER 



"Data Set Number - " " , DDD , 2X , 1 4A " ; J , To 1 d$ 

" T v 1 Tv2 Tv3 Tldl T 1 d2 TldZ 
PRINT USING " 10X ,8( MOD. DO, 2 X )“ ; T< 0 >,T( 1 ) , T< 2 ) , T< 3 ) , T( 4 ) , T( 5 ) , Tvav , Tlev 



3096 


PRINTEP IS ' 


701 


3100 


PRINT 




3104 


PRINT USING 


“10X 


3108 


PFINT 




3112 


PRINT USING 


“10X 


Tvav 


Tldav 




3116 


PRINT USING 


"10X 


3120 


PRINT 




3124 


PRINT USING 


" 6 X , 


H 


The tat " 




3128 


PRINT USING 


" 6X , 


2 ) ( L 1 / m 


" 2 . K > < K > " ■ 




3129 


IF lpos>0 THEN 


3130 


J j*0 




3132 


FOR Ji-0 


TO 5 






Ulall Temperatures (Deg C) Tnave Odp 

2 3 4 5 6 Dec C > (U/m* 
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3133 TpiJi**Tvlpos*E-*Jj+6) 

3134 Jj-Jj+1 

3136 NEXT J; 

3137 T nn = 1 pos+ 1 

3138 PRINT USING "6X ,D, IX ,7( MDO.DO ) , IX .2^ M2 .3DE ) , 1> , 1 < MOD. DC 1 . Tnn. Tp C 
, Tp ( 1 ) , Tp ( 2 ),Tp 3 ) , Tp 4 },Tp(5/,Twe' Ipos ) , Odp< Ipos , HtuDe ( I pot , The t ab Ipos ) 

3139 GOTC 3177 

3140 END IF 

3141 J j - 0 

3143 FDR Ii«0 TO N h t * 1 

3143 F OF Ji-0 TO 5 

3144 Tp ( J i )«T { 1 1 • 5+ J j + 6 ) 

3146 Jj«Jj+l 

3152 NEXT Ji 

3 15B Jj=li+1 

3160 F OF J i-e TO 4 

3164 Tn( J : >*1+Ji 

31 BB NEXT Ji 

3172 PRINT USING “ 6X , D , 1 X , 7( MDD. DO ) , 1 X , 2 ( MZ . 3DE ) , 1 X , 1 < MDD . DD ) “ . Tn( I : ) , T 

p( 0 ),Tp( 1 ) ,Tp( 2 ),Tp( 3 ),Tp(4 ),Tp(5 ) ,Tua( 1 1 ) ,Gdp< 1 1 ),Htube< 1 1 >,Thetab( Ii ) 



3176 NEXT Ii 

3177 PRINT 

3179 PFINT USIN6 *6X,*“ Heat Flux and Tdel fron Church: 1 1 /Chu Correlation i 
s “M(MZ.3DE ),2X,1(MDD.DD) - ;0ch, The tab ( Ipos > 

3180 PRINT 

3182 Dk * 1 

3184 IF In-0 THEN 

3188 BEEP 

3192 INPUT “OK TO STORE THIS DATA SET ( 1 - Y( de f au 1 1 ) , 0-N )? * , Dk 

3196 END IF 

32001 J«the counter for data sets 

3204 IF Dk* 1 OR In-1 THEN J*J+ 1 

3208 IF Ok * 1 AND In-0 THEN OUTPUT &F 1 1 e 1 ; Enf ( « ) , Vol t ( • ) , Anp ( * ) 

3212 IF In-1 DR Dk - 1 THEN OUTPUT ©Plot . Qdp( * ) ,Htube( * ) , Thetab( * ) 

3216 Go_on«l 

3220 I F~ I n*0 THEN 

3224 BEEP 

3228 INPUT "UILL THERE BE ANOTHER RUN ( 1 «Y ( de f au 1 1 0-N P", Gc_on 

3232 Nrun«J 



3236 IF 6o_on-0 THEN 3272 

3240 IF Gc_on<>0 THEN Repeat 

3244 ELSE 

3248 IF J<Nrun+l THEN Repeat 

3252 END IF 

3256 St-1 

3260 BEEP 

3264 IN D UT “ARE YOU SURE YOUF READY TO TERMINATE < 1 »Y ( DEFAULT ), 0-NO P ", St 

3268 Gc_on« 1 

3272 IF St>e THEN 3280 

3276 IF St-0 THEN GOTO 3240 

3280 IF In-0 THEN 

3284 BEEP 

3288 PRINT 

3292 PRINT USING “ 1 0X , " “NOTE "" ,32, M " data runs were stored ;r file “ ’ 

, 10A“ ; J-l ,D2f il-5 

3296 ASSIGN ©Filel TO * 

3200 OJTPU T @F i le2 ; N-ur- 1 

3304 ASSIGN ©Filel TC DlfileS 

33CS ENTER l 1 el ; Da t eS , Ldt c < * ) , 1 1 1 , Bop , Nnt , Ne tp . Nr * , Co r ' , ‘ 

3312 OUTPUT @F 1 1 e2 ; Det eS , Ldt c ( * ) , 1 1 1 . Bop , Nht , No tp , Nr t , Cor •- 
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FOF 1*1 TO N-ur-i 



ENTEF : 1 e 1 . E ^ f ! • > , Uc * t • * i . Amp • • 



OUTPUT @F : i e :,c n f ( . , , L- c 1 1 < * /.hmp' * 



NEXT I 



ASSIGN ©Filel TG • 



j - 36 
3340 
3344 

334e 

3352 

i 1 e 

3356 

2360 

3364 



PURGE ■OUMnY" 

END IF 

BEEF 

PRINT 

PRINT USIN6 ' 1 0X , " “ NOTE " " ,Z2 ” X-Y pairs were sto-ed in plot oata f 



, 1 0A “ i J- 1 , P f : 1 eS 



ASSIGN @F i 1 e2 TO • 
ASSIGN ©Plot TO * 
BEEP 



3368 SUBENO 
3372' 

3376 IXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
2380 1 

3384 DEF FNKcu(Tcu ) 

3388' OFHC COPPER 

3382 Tk -Tcu+273 . 1 5 1C TO K 

3396 Kcu*434- . 1 1 2 * Tk '250-300K USE FOR R-l 14 ©2.2 C 

34001 Kcu-433.0-. 1 * T k ' 200-400K USE FOR R-l 13 @47.5 C 

3404 RETURN Kcu 
3408 FNENO 
3412' 

3416 OEF FNhu< T ) 

3420! CURUE FIT OF UISCOSITY 
3424 Tk-T+273.15 'C TO K 

3426 Mu*EXP< -4 . 4636+< 1011 .47/Tk > )♦ 1 . 0E-3 'R-114 170-360 K 

3432! hu«. 00001 34 «( 1 0^(503/ (Tk-2. 15 > > > «R1 13 

3436 RETURN Mu 

3440 FNENO 

3444 I 

3448 DEF FNCp(T) 

3452' CURUE FIT OF Cp 

3456 Tk *T + 272 . 15 'C TO Y 

3460 Cp*. 40188+1 .6S007E-3»Tk*l . S1494E-S*Tk *2-E .67853E-10*Tk'3 > R- 1 1 4 180-400 K 

34641 Cp*( S25 4 - 1 .C3*T )• . 001 IR-113 

3468 Cp*Cp*1000 

3472 RETURN Cp 

3476 FNENO 

3480' 

3484 DEF FNRho< T > 

3488 Tk«T^273.15 'C TO K 

3492 X«1-(1.8*U/7S3.95> 'K TC R 

34 96 Ro'36. 32 + 61 . 1 464 1 4 *x * ( i /3 )+ 1 6 . 4 1 801 5* X + 1 7 . 476858* X " . 5+ 1 . 1 1 9828* X* 2 
3500 Ro'Ro/ . 0E2426 'R-114 

2504 i Ro* 1 .6207479E + 3-T*(2.2l8634E + T*2.3578291E-3 > ' F - 1 1 3 

2508 RETURN Ro 
3512 f n c n0 
2516' 

2520 DEF FNPr(T) 'GOOD FOR R-l 14/R-l 13 
3524 Rr*FNCp( T )*FNMu r T >/FNK< T > 

2528 RETURN Pr 
3532 FNENO 



3536' 

3540 OEF FNK ( T ) 

3544 i T <360 K WITH T IN C 
2549 K* . 071 - . 00026 1 * T 
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3552 PET URN y 
1 5 5c FNEND 
3560' 

3564 DEF FNTanh(Fe ) 

3560 P-E>P(Fe) 

3573 O-EXP(-Fe) 

3576 Tanh«< F-0 >/<P+Q ) 

3580 RETURN Tanh 
3584 FNEND 
3588' 

3593 OEF FNTvsv(U) 

3596 COM /Cc/ C( 7 > 

3600 T-C(0) 

3604 FOR >1 TO 7 
3608 T-T + C(I)«VM 
3612 NEXT I 
3616 RETURN T 
3620 FNENO 
3624 1 

3628 DEF FNBeta ( T ) 

5632 Rop-FNRho( T+ . 1 ) 

3636 Ror>-FNRho ( T- . 1 ) 

3640 Beta--2/(Rop+Ron )• ( Rop-Rom )/.2 

3644 RETURN Beta 

3648 FNENO 

3652 DEF FNPoly<X ) 

3656 COM /Cply/ A< 10 , 10 ) , C( 10 ) , B( 4 ) , Nop t Iprnt ,0po , I log 
3660 Xl-X 
3664 Foly-B<0) 

3668 FOR 1-1 TO Nop 

3672 IF IIoq-1 THEN X 1 »L06( X ) 

3676 Poly-Poly+B< I )»XKI 

3680 NEXT I 

3684 IF IIoq- 1 THEN Po 1 y-EXP ( Po 1 y ) 

3688 RETURN Poly 
3692 FNEND 
3696' 

3700 IXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
3704 I 

3708 SUE Poly(Ofile$(* ) , Np , 1 1 n ) 

3712 Dir R( 10 ),S( 10 >,Sy< 12 ) ,S*< 12 ) ,X*< 100 ),Yy( 100 ), Xy< 17 ) 

3716 COM /Cply/ A< 10 , 10 ) , C( 10 ) , B< 4 ) , N , Iprnt , Opo, I log 
3720 COM /Xxyy/ Xp(5),Vp(5) 

3724 FOR 1-0 TO 4 
3728 B(I>«0 

3732 NEXT I 
3736 In-1 
3740 BEEF 

3744 INPUT "ENTER OAT A FILE NAME " , D f 1 1 e$ ( 0 ) 

3746 BEEP 

3752 INPUT "ENTEF NUMBER OF X-Y PAIRS", Np 
3756 BEEP 

3760 INPUT “LIKE TO EXCLUOE OATA PAIRS (1-Y,0-N DEFAULT ) ) 9 " , 1 ed 
3764 IF Ied-1 THEN 
3768 BEEP 

3772 I N d UT "ENTEF NUM5EF OF PAIRS TC EE EXCLUDED" , Ipe> 

3776 END IF 

3780 ASSIGN §>F:lc TO DfileS(0> 

3784 N = 2 
3786 BEEP 



163 



3792 

3' 7 96 

3600 

3904 

3806 

3612 

3616 

3620 

3624 

3626 

5632 

3836 

3840 

3844 

3846 

3652 

3656 

3860 

3864 

3868 

3872 

3876 

3880 

3884 

3888 

3892 

3896 

3900 

3904 

3908 

3912 

3918 

3320 

3924 

3928 

393 ^ 

3536 

3940 

3944 

3946 

3952 

3956 

3960 

3964 

3966 

3572 

3976 

3980 

3984 

3988 

3992 

3996 

4000 

4004 

4006 

4012 

4016 

4020 

4024 

4026 



INPUT "ENTEF THE OROEF OF POLYNOMIAL (0EFAULT*2 *' , N 
FOF 1-0 TO N 
Sy< 1 )*0 
Sx< i >«e 
NEXT I 

IP I ed* 1 ANO Im«l THEN 
FOR 1*1 TO Ipe> 

ENTER ©File.Xyf * ) 

NEXT I 
ENO IF 

FOR 1*1 TO Np-Ipex 

ENTER @F i le ; Xy( * ) 

IF Opo*0 THEN 
Y«Xy( 1 1 n- 1 ) 

X*Xy ( 11+Itn) 

ENO IF 

IF Opo* 1 THEN 
Y*Xy( 5+Itn ) 

X«Xy( 11+Itn ) 

ENO IF 

IF Opo«2 THEN 
Y*Xy( 5+Itn ) 

X *X y ( Itn-1 ) 

END IF 

IF I log*l THEN 
X*L06( X ) 

Y*L06( Y ) 

ENO IF 
Xx( I )*X 
Y y ( I )«Y 
R ( 0 )*Y 

Sy ( 0 )*Sy ( 0 ) + Y 
5 ( 1 )«X 

Sx( 1 )-Sx( 1 ) + X 
FOR J-l TO N 

R( J >*R< J-l )*X 
5y( J )*Sy ( J ) + P< J ) 

NEXT J 

FOR J*2 TO N*2 
S' J )*S( J-l )*X 

Sx( J )«Sx( J ) + S( J ) 

NEXT J 
NEXT I 
Sx ( 0 )*Np 
FOR 1*0 TO N 
C( I )*Sy( I ) 

FOR J-0 TO N 

A( I , J >*S-» ( I + J ) 

NEXT J 
NEXT I 

FOR 1*0 TO N-l 

CALL Di vi de< I ) 

CALL Subtrac t ( I * 1 ) 

NEXT I 

B ( N )»C( N )/A( N,N ) 

FOR 1*0 TO N-l 

B( N-l-I ;*C< N- 1 - 1 ; 

FOP J«0 TO I 

B( N-l -I i=B( N-I-I ‘ - N- 1 - 1 , N- J ) • B f N- J ) 

NEXT J 
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4032 B( N- ] - : »-B( N-l-I A i N— 1 — I , M- ] - I 

4036 NEXT I 

4040 IF Iprnt*C THEN 

4044 PRINT USING "12> . " "EXPONENT COEFF I C I ENT " " ” 

404 8 FQF >C TO N 

405: PFIN t USING " 1 5X , DC , 5X , MO . 70E “ . I , B( I ) 

4056 NEXT I 

4050 PRINT ” " 

4064 PRINT USING ’ 12X , " ’DATA POINT > Y Y< CALCULATED < 

SCREPANCY*'*" 

4066 FOR 1-1 TO Np 

4072 Yc-BiO) 

4076 EOF J«I TO N 

40B0 Yc»Yc+B( J ) • X >: ( I J 



40B4 

4066 

4092 

4096 

4100 

4104 

4108 

4112* 

4116 

4120 

4124 

4 1 2 B 

4132 

4136 

4140 

4144 

4146 

4152 

4156' 

4160 

4164 

4 1 6B 

4172 

4176 

41B0 

4184 

4186 

4192 

41961 

4200 

4204 

4206 

4212 

4216 

4220 

4224 

4228 

4232 

4236 

4240 

4244 

4246 

4252 

4256 

4260 

4 264 



NEXT J 
D'Yy( I )-Yc 

PRINT USING "15X,30,4X ,4( M0.50E, IX >■ , I t Xx( I >,Yy( I\Yc,D 
NEXT I 
END IF 

ASSI6N GFiIe TO ♦ 

SUBENO 

SUE Oivide(M) 

COn /Cply/ AJ 10, 10 ),C( 10 ) ,B(4 ),N, Iprni ,0po, I log 
FOR I*M TO N 
Ao*=A( I , r ) 

FOR J-r TO N 

A( I , J )«A( I , J )/Ao 
NEXT J 

C( I >*=C< I )/Ao 
NEXT 1 
SUBENO 

SUE Subtract ( K ) 

COM / Cply/ A( 1 0 , 1 0 ) t C ( 10), 6(4 ) , N , I prnt , Opo , 1 1 og 
FOR I C K TO N 

FOR J ■ K - 1 TO N 

A( I , J >«A( K - 1 , J >-A( I , J ) 

NEXT J 

C( I >*C( ► -1 >-C( I ) 

NEXT I 
SUBENO 

SUB PI in 

COM /Cply/ A( 10, 10 ) ,C( 10 ),E( 4 ),N, Iprnt ,0po, I log 
COM /Xx yy/ X*(5),Yy<5) 

PRINTER IS 705 
BEEP 

INPUT ’SELECT <0«h/h02 sane tube , l*h( HF )/h( sn ) " , Irt 
BEEF 

INPUT " UHI CH Test ( 1-6 . 7 , 0=-2 . 2 ) M sat 
Xnn B 0 
Xnax c 1 0 
Xst ep = 2 
IF I r t *0 THEN 
Yni n-0 
Ynax'l . 4 

V 5 tep= . 2 

ELSE 

Y r i r = 0 



01 
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4266 

4272 

427E 

4280 

4284 

4288 

4292 

4296 

4 300 

4 304 

4306 

4312 

4316 

4320 

4324 

4326 

4532 

4336 

4340 

4344 

4348 

4352 

4356 

4360 

4364 

4368 

4372 

4376 

4380 

4384 

4388 

4392 

4396 

4400 

4404 

4408 

4412 

44 IE 

4420 

4424 

4428 

4432 

4436 

4440 

4444 

4448 

4452 

4456 

4460 

4464 

4466 

4472 

4476 

448C 

4484 

4468 

4452 

44 56 

4S0O 

4S04 



Yma * m IS 
Y 5 tep*S 
END IF 
6EEF 

PRINT " IN ; S P 1 . IF 2300,2200,6300,6800." 

PRINT "SC 0, 100, e, 100 . TL 2,0.' 

Sf>*100/( Xma>-Xmin ) 

S f y * 1 00 / ( Y ma > - Y ru n ) 

PRINT "PU 0,0 PD" 

FOR Xa»Xr»un TO Xna> STEP Xstep 
X*( Xe-Xmi n )»Sf * 

PRINT “PA*\X, “ ,0, XT. " 

NEXT Xa 

PRINT "PA 100,0. PU. M 
PRINT ‘PU PA 0,0 PD** 

FOP Ya*Ymin TO Yna- STE P Ystep 
Y-< Ya-Yni n XSf y 
PRINT “ PA 0, " ; Y, "YT" 

NEXT Ya 

PRINT ‘PA 0,100 TL 0 2 ” 

FOR Xa*Xnin TO Xma> STEP Xstep 
X» ( Xa-Xm n >* S f x 
PRINT "PA" i X , 100. XT" 

NEXT Xa 

PRINT "PA 100,100 PU PA 100,0 PD" 

FOR Ya-Ynin TO Ynax STEP Ystep 
Y»( Ya-Ynin )*Sf y 
PRINT ‘PD PA 100,",Y,"YT" 

NEXT Ya 

PRINT ‘PA 100,100 PU" 

PRINT ‘PA 0,-2 SR 1.5,2" 

FOR Xa a Xmin TO Xmax STEP Xstep 
X*( Xa-Xnin XSf> 

PRINT "PA" ; X , " , 0 . 

PRINT "CP -2,-l.LB‘iXa. "" 

NEXT Xa 

PRINT "PU PA e,0" 

FOP Ya*Yr. in TO Yma a STEP Ystep 
IF AES ( Ya XI . E- 5 THEN Ya-0 
Y*(Ye-Yrun)*Sfy 
PRINT ‘PA 0 , ‘ , Y , “ 

PRINT "CP -4 , - . 25 ; LB " ; Ya . 

NEXT Ya 

XlabelS*"Oil Percent “ 

IF Irt-0 THEN 

Ylabel5*“h/h0X" 

ELSE 

YlabelS="h/hsmooth" 

END IF 

PRINT "SF 1 . 5 , 2 . PU PA 50,-10 CP *' . -LEN( X label S )/2 . ‘ C . LB ** ; X 1 abe 1 S . " “ 
PRINT "PA -11, SC CF 0, " ; - LEN( Y labels )/2*S/6 ; "DI 0 , 1 ; LB " ; Y 1 abe 1 $ , " ‘ 
PRINT "CP 0,0" 

I p n«C 
EEE P 

INPUT "WANT TO PLOT DATA FROM A PILE ( 1 , 0«N )? " , 0* p 
I C n = C 

IF 0 \ p-l THEN 
BEEP 

INPUT "ENTER THE NAME of the DATA FILE" ,D_f l ieS 
BEEP 
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45C6 

45:2 

*151 E 

4 52C 

45:-: 

4526 

4532 

4536 

454C 

4544 

4548 

455: 

4556 

4560 

4564 

4568 

4572 

4576 

4580 

4584 

4588 

4592 

4596 

4600 

4604 

4608 

4612 

4616 

4620 

4624 

4628 

4632 

4656 

4640 

4644 

4648 

4652 

4656 

4660 

4664 

4668 

4672 

4676 

4680 

4684 

4686 

46 52 

4696 

4700 

4704 

4708 

4712 

4716 

4~20 

4724 

4726 

4722 

4"3S 

4 74 C 

4744 



INPUT "SE^EC" < 0 = L I NEAF. , l-.06( X 1 lo; 

ASSIGN ©File TC C_'i leS 
6EEF 

INPUT "ENTEP THE BEGINNING RUN NUMBER " , Md 
EEEF 

INPUT "ENTEF THE NUM6EF OF X-Y PAIRS STOREO " , Npa 1 r s 
BEEP 

I NPUT "ENTER OESIRED HEAT FLUX",G 
BEEP 



" 1 
"5 
“5 
"6 
“7 



PRINTER 16 1 
PRINT USING " 4 X , 1 
PRINT USING " 4 X , ' 

PRINT USING " 4 X , ' 

PRINT USING " 4 X t ' 

PRINT USING “4X , ' 

PRINT USING ’ 4X , ' 

INPUT Sym 

PRINTER 15 705 

PRINT "PU or 

IF Syr** 1 THEN PRINT ’SM*’ 

IF Syn-2 THEN PRINT "SM+" 

IF Syw-3 THEN PRINT ’SMo’ 

Nn-4 



Select e symbol " " 



Star 2 Plus s ign* " " 
Circle 4 SQua^e’"" 

Rombus * " “ 

Right-side-up triangle""" 
Up-s 1 de-down triangle"’" 



IF Ilog*l THEN Nn* 1 
IF Md> 1 THEN 

FOR 1-1 TO (hd-1 ) 

ENTEP @File.Xa,Ya 
NEXT I 
END IF 
01-0 

IF Ilog-1 THEN O-LOG(O) 
FOR 1*1 TO Npairs 

ENTER ®F 1 1 e . Xa , 8 ( * ) 

Y a*B ( 0 ) 

FOR K* 1 TO Nn 



Ya-Ya + 8 ( K >*Q*K 
NEXT K 

IF 1 ]og*l THEN Ye-EXP( Ya ) 

IF Ilog-0 THEN Ya-Ql/Ye 
IF Irt-0 THEN 
IF Xa-0 THEN 
Yo-Ya 
Ya* 1 
ELSE 

Ya*Ya/Yo 
END IF 
EL.SE 

Hsm=FNHsmocth(C,Xa,lsat ) 
Ya=Ya/Hsm 
END IF 
X*( 1-1 >«Xa 



Yy( 1-1 )- Ya 
X«< Xa-Xmin )*Sf x 



V * ' Y e - V r 1 n > * S f y 
IF Sym/3 THEN PRINT 
IF Symv4 THEN PRINT 
PRINT "PA" ,Y , Y, 

IF Syr 3 THEN PRINT 
1 P Syr* 4 THEN PRINT 
IF S y r = 5 THEN PRINT 



"SM" 

"SF 1. 4,2.4 " 

"SR 1.2, 1.6“ 
" UC2 .4.99,0, 
"UC3 , 0 , 99 , -3 



8, -4,0,0. 6,4, C 
-E. -3, 6, 3, 6,3. 
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4' , 48 

4?b: 

4756 

4760 

4764 

4766 

4777 

4776 

4780 

4784 

4788 

4797 

4756 

4800 

4804 

4808 

4612 

4ei6 

4820 

4824 

4828 

4632 

4836 

4840 

4844 

4848 

4852 

4856 

4860 

4864 

4868 

4672 

4876 

4880 

4884 

4888 

4 652 

4895 

4900 

4904 

4906 

4912 

4916 

4920 

4924 

4928 

4932 

4936 

4940 

4944 

4948 

4952 

4956 

4960 

4964 

496e 

4972 

4976 

4980 

4964 



IP Sy n-E THEN PRINT " uCC , 5 . 3 , 99 , 3 , -8 , -6 . C , 3 , 8 . 

IP Sy" a ~ THEN PRINT "UCC. -5. 2, 99, -3, 6, 6,0. -3, -6. " 

NEXT I 

EEE C ' 

ASSI6N le TO • 

END IF 

PRINT "PU sir* 

B6EP 

INPUT "WANT TO PLOT A PDLYNDMIAL U «Y , 0-N ) 7 M , D1 p 
IF 0>p«l THEN 
BEEP 

INPUT "SELECT < 0-L INEAR , 1 «LDG< X , Y ) ) “ , 1 1 og 
Iprnt - 1 

CALL Poly( Itn ) 

FOP Xe-Xnin TD Xma* STE P Xstep/25 
Icn« Icn + 1 
Ye-FN p ol y ( Xa ) 

Y-( Ye-Ynin )*S f y 
X = ( Xa-Xmin )*$f > 

IF Y <0 THEN Y-0 

IF Y> 100 THEN GOTD 4868 

Pu-0 

IF lpn-1 THEN Idf-Icn HDD 2 
IF Ipn-Z THEN Idf-Icn MDD 4 
IF Ipn-3 THEN Idf-Icn MOD 8 

IF Ipn-4 THEN Idf-Icn MOD 16 

IF Ipn-5 THEN Idf-Icn MDD 32 

IF Idf-1 THEN Fu-1 

IF Pu-0 THEN PRINT * PA “ , X , Y , * PD " 

IF Pu-1 THEN PRINT " PA * , X , Y , “ PU " 

NEXT Xa 
PRINT "PU * 

I on* I pn+ 1 
60TD 4480 
END IF 
BEEP 

INPUT "WANT TO QUIT * 1 = Y , 0-N ) 7 “ , I qu 1 1 
IF Ipuit-1 THEN 4904 
GOTO 4480 
PRINT "PU SP0“ 

SUBEND 

SUB Stats 

PRINTER IS 701 

J-0 

K-e 

BEEP 

INPUT "PLDT FILE TO ANALYZE 7 " , F i ] eS 

ASSIGN ©rile TO Files 

BEEP 

INPUT "LAST RUN No 7 ( 0 = QUI T ) " , Nn 

IF Nn=0 THEN 5092 

N n - N n - J 

Sr-0 

Sy-0 

Sr-0 

Sxs«0 

5 y £ = 0 
Sz £-0 

p DP 1=1 TO Nn 
J-J + l 
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493S ENTEF @R 1 1 e . C , T 
49°: H*0 '1 

4 99E 5..*S^Q 

5 000 S^-S*s*Q*2 

5004 Sy-5y*T 
500B Sy i- Sy s + T 2 
5012 Sz-Sz+H 
5016 Sz s-Sz s + 2 

5020 NEXT 1 
5024 Gave-S/./Nn 
502B Tave-Sy/Nn 
5032 Have»Sz/Nn 

5036 5devq-5QP( AB5( < Nn* 5 a 5-S* " 2 >/<Nn«<Nn-l ) ) ) 

5040 Sdevt«SQR( ABS< <Nn«Sys-S> "2 )/<Nn*( Nn-1 ) > ) ) 

504 4 Sdevh-S0R< ABS< ( Nn*Szs-Sr "2 )/< Nn*( Nn-1 ) ) ) ) 

5040 Sh- 1 00* Sde vh/Heve 
5052 Sq- 100* Sde vq/Qave 
5056 5 t - 1 00* Sde vt /Tave 

5060 IF K-l THEN 50B4 
5064 PRINT 

506B PRINT USING “11X, ““DATA FILE UA“ ; F i let 
5072 PRINT 

5076 PRINT USING MIX, ""RUN Htube SdevH Odp SdevQ Thetat SdevT““ 

5080 K-l 

50B4 PRINT USING * 1 1 X , DO , 2< 2X , D . 3DE , 1 X , 3D . 2D ), 2X , DO . 3D , 1 X , 3D . 20' ; J , Have , Sh , Qave 

, Sq , Ta ve , St 

50B8 GDTO 4840 

5092 ASSIGN eFilel TO * 

5096 PRINTER IS 1 
5100 SUBEND 
5104 SUB Coe f 

510B COM /Cply/ A( 10 , 10 ),C( 10 ) ,B<4 ),N, lprnt ,0po, I log 
5112 BEEP 

5116 INPUT "GIUE A NAME FOR CROSS-PLDT FILE“,Cpf* 

5120 BEEP 

5124 INPUT "OUTPUT TYPE <0«q vt Ot , 1-h vs Ot , 2-h vt q)*’,0pc 
512c CREATE BOAT CpfS,6 
5132 ASSI6N ©File TO CpfS 
5136 BEEP 

5140 INPUT “SELECT ( 0-LI NEAR , 1 »LD6( X , Y ) >“ , 1 log 
5144 BEEP 

514B INPUT “ENTER OIL PERCENT ( - 1 -STDP ) “ , Bop 
5152 BEEP 

5 15E INPUT “ENTER TUBE NUMBER (1, 2, 3, 4, OR 5>‘\Itn 
5160 IF 6op 0 THEN 5176 
5164 CALL PolyC Itn) 

5166 OUTPUT @F i 1 e . Bop , B ( * ) 

5172 GOTO 5144 

5176 ASSIGN ©Pile TC • 

5180 SUBEND 
51B4 i 

51 B8 ' XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

5192' 

5 1 SB SUE Plot 

5200 cor /Cply/ A:lC,!C , C 1 0 1 , E ( 4 ) , Nop , 1 prnt , Opo , 1 1 eg 
5204 oin x y ; 17 ) 

5206 INTEGER U 
5212 PRINTER IS 1 
5216 BEEP 
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s::c 

5224 

5:20 

5:32 

5236 

5240 

5244 

5248 

5252 

£256 

5260 

5264 

5266 

5272 

5276 

52B0 

5264 

52B6 

5292 

5296 

5300 

5304 

5306 

5312 

5316 

5320 

5324 

5326 

5332 

5336 

5340 

5344 

5346 

5352 

5356 

5360 

5364 

5368 

5372 

5376 

5380 

5364 

5368 

5392 

5396 

5400 

5404 

5408 

5412 

5416 

5420 

5424 

542e 

5432 

5436 

5440 

5444 

5446 

545: 

5456 



I dv* 1 

INPUT "LIKE DEFAULT VALUES FOF PLOT il-X DEFAULT ,0«N n ' ,12, 

Opo-C 

BEEP 

PRINT U5IN6 ’4X, "Select Option "" 

PRINT USING ’ 6 X , " " 0 0 versus del ta-T< DEFAULT ) 

PRINT USING ’ 6 X , * “ 1 h versus delta-T""" 

PRINT USING ‘6X , ’ “ 2 h versus q " " " 

INPUT Opo 
6EEP 

INPUT “SELECT UNITS ( 0-SI ( DEFAULT ), 1 -ENGLI SH I un 
PRINTER IS 705 
IF IdvOl THEN 
BEEP 

INPUT ‘ENTEF NUMBEP OF CYCLES FOR X-AXIS" ,Ca 
BEEP 

INPUT 'ENTER NUMBER OF CYCLES FOR >-AXIS H ,Cy 
BEEP 

INPUT "ENTER MIN X-VALUE (MULTIPLE OF 10)",Xmn 
BEEP 

INPUT 'ENTER MIN Y-VALUE (MULTIPLE OF 10>',Ymn 
ELSE 

IF Opo-0 THEN 
Cy-2 
C**2 
Xmin* 1 
Ynm*1000 
END IF 

IF Opo* 1 THEN 
Cy*2 
Cx = 2 
Xn in* 1 
Ymin-100 
ENC IF 

IF 0po=2 THEN 
Cy-2 
C> *2 

Xfi 1 n= 1 000 
Yni n= 100 
END IF 
END IF 
BEEF 

PRINT "IN;SP1;IP 2300,2200,6300,6800.' 

PRINT 'SC 0,100,0,100.10 2,0; ' 

Sf>-100/C* 

S f v = 1 00/Cy 
EEEP 

INPUT “WANT TO BY-PASS CAGE (1*Y, 0*NO( DEFAULT >“, I bp 
IF I bp * 1 THEN 5906 
PRINT "PU 0,0 PD" 

Nn*9 

FOR 1=1 TO Ca+1 

Xat*Xnin«10"( I - 1 ) 

I R I*Cx+I THEN Nn= 1 
POP J= 1 TO Nn 

IF J«I THEN PRINT - Ti_ : ?“ 

IF J=2 THEN D FIN T “Tl 1 0' 

> e * X a t * J 

X =LGT( Xa/Xni n )* S f >■ 

PRINT “ P A C , V T . “ 
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5«i6C 

5464 

SAGS 

547: 

5*176 

546C 

5484 

54S& 

54S: 

54 96 

550(2 

5504 

5508 

5512 

5516 

5520 

5524 

5526 

5532 

5536 

5540 

5544 

5548 

5552 

5556 

5560 

5564 

5568 

5572 

5576 

5580 

5584 

5588 

5592 

5596 

5600 

5604 

5608 

5612 

5816 

5620 

5624 

5628 

5632 

5636 

5640 

5644 

S64e 

5852 

5656 

5860 

5664 

5568 

5672 

5678 

5680 

55 64 

sees 

5852 

5896 



NEXT J 
NEXT I 

PRINT "PA 10C ,C . PL 1 
PRINT "PU PA C,0 PO" 

Nn«S 

c OF l«I TO Cy+1 

Yat*Yrun« 10" ( I-l ) 

I r I«Cy+l THEN Nn-1 
POP J=1 TO Nn 

IF J-l THEN PR I N T " TL 2 0 " 

IF J *2 THEN PRINT " TL 1 0 " 

Y«LGT(Ya/>nin)*5fy 
PRINT "PA 0, • ; Y , “YT“ 

NEXT J 
NEXT I 

PRINT "PA 0,100 TL 0 2" 

Nn-9 

FOR 1-1 TO C* + l 

Xat-Xrun* 10" ( 1-1 ) 

IF I-Cx+1 THEN Nn-1 
FOR J-l TO Nn 

IF J-l THEN PRINT " TL 0 2" 

IF J>1 THEN PRINT "TL 0 1” 

Xa-Xat • J 

X-L6T<Xa/Xr*un )*Sf x 
PRINT " PA * ; X , " , 100 ; XT- 
NEXT J 
NEXT I 

PRINT "PA 100,100 PU PA 100,0 PD" 

Nn* 9 

FOR 1-1 TO Cy+1 

Yat-Ynin«10"< 1-1 ) 

IF I-Cy+1 THEN Nn*l 
FOP J-l TO Nn 

IF J-l THEN PRINT “TL 0 2" 

IF y, 1 THEN PRINT M TL 0 1 " 

Ya-Yat*J 

Y - L G T ( Ya/Ynin ; • 5 f y 
PRINT "PD PA 1 00 , " , Y , " YT " 

NEXT J 
NEXT I 

PRINT "PA 100,100 PU" 

PRINT “PA 0,-2 SR 1.5,2" 

I i -LGT < Xpu n ) 

F OP 1-1 TO Cx+1 

Xa-Xmn*10"( 1-1 ) 

X-LGT ( Xa/Xmn )*Sfx 
PRINT " P A " ; X , " , 0 . " 

IF I i >-C THEN PRINT "C^ -2 , -2 ; LB 1 0 . PP -2 , 2 ; LB " ; 1 1 ; " “ 
IF 1 1 <0 THEN PRINT "CP -2 , -2 ; LB 1 0 ; PR 0 , 2 ; LB “ ; 1 1 ; " “ 
Ii«Ii+l 

NEXT I 

PRINT "PU PA C,0" 

Ii-LGT(Yrnin) 

Y 1 0 * 1 0 

RO c 1-1 TO Cv+l 

Ya- Ym n* 1C " ( 1-1 ) 

Y*LGT( Y a / Y n i n ; ♦ 5 f y 
PRINT "PA 0 , " ; Y , " ” 
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5700 PFINT -CP -4,-.25.^ElC.FF -1 , 2 . LB " . i i . “ " 

5704 Ii-Ii+i 

s~e6 next i 
5^:: beep 

5 *'16 Idl-1 

572C INPUT ‘UANT USE DEFAULT LABELS ( 1 -Y( DEFAULT ), 0-N >*> " , I d 1 
5724 IF IdK'i THEN 
5^26 BEEP 

5-32 INPUT "ENTER Y-LAEEw” , XiabeiS 
5"36 BEEP 

5740 INPUT "ENTER Y-LABEL" , V 1 abe IS 

5744 END IF 

5746 IF Ddo<2 THEN 

5752 PRINT ”5R 1,2. PU PA 40,-14;" 

5756 PRINT "LB( T . PF. -1.6,2 PD PR 1.2,0 PU ; PF . 5 , -4 , LBuo . PP .5,1;" 

5760 PRINT "LE-T ; PR . 5 , - 1 ; LE so t ; PR .5,1." 

5764 IF Iun-0 THEN 

576B PRINT "LE > / ( K ) ‘ 

5772 ELSE 

5776 PRINT "LB ) / < F ) " 

5780 END IF 

5794 END IF 
5788 IF Opo-2 THEN 
5732 IF Iun-0 THEN 

5796 PRINT "SR 1.5,2;PU PA 40,-14;LBq / (U/n;SR 1 , 1 . 5 , PR 0.5,1; LB2 ; SR 1 

. 5 , 2 ; PR 0. 5,-1. LB)" 

5B00 ELSE 

5804 PRINT "SR 1 . 5 , 2 ; PU PA 34,-14;LBq / (Btu/hr;PR .S,.S;LB.;PP .5, -.5; 

5808 PRINT "LB f t ; PR .5,1:SR 1,1.5; LB2 ; SR 1.5, 2. PR .S,-1;LB>." 

5812 END IF 

5816 END IF 
5B20 IF Opo-0 THEN 
5824 IF Iun-0 THEN 

5826 PRINT "SP 1 . 5 , 2 ; PU PA -12,40.DI 0,I;LBq / (U/n.PR -1,0. 5. SR 1,1. 5 ; L 

B2.SR 1 . 5 , 2 ; PR 1,.5;LB>* 

5S32 ELSE 

5636 PRINT "SR i.5,2;PU PA -12,32. D2 0,1. LBq / (etu/hr,PR - . £ , . 5 ; LB . . PR 

.5, .5. " 

5840 PRINT "LBft.SR 1,1. 5. PR - 1 , . 5 ; LB2 ; PR 1,.5.SR 1.5,2;LB)“ 

5644 END IF 

5848 END IF 
5852 IF Opo>0 THEN 
5656 IF Iun-0 THEN 

5860 PRINT "SP 1 . 5 , 2 ; P U PA -I2,38;DI 0,1. LBh / OJ/n.PF -1,.5 SF 1,1. 5. LB 

2.SP 1 .5,2. PR .5, .5; " 

5864 PRINT “LE. .PR .5,0. LEp )" 

sese ELSE 

5872 PRINT "SP 1.5,2;PU PA -12,28;DI 0,1. LBh / (Btu/hr.PR - . 5 , . 5 . LE . ; PF 

. 5 , . 5 ; " 

5E76 PRINT "LBft.PR -1..5.SR 1,1.5.LE2.SR 1.5,2;PF .5,.5;LB..PR .5, .5; 

L B F ) ” 

5880 END IF 

5864 END IF 
5688 IF Idl-e THEN 

5892 PP I NT "SF 1 . E , 2 ■ PU FA 5C.-16 CP" . -LEN> > labels i/2 . "0 . LB " ; X I abe 1 $ . " " 

5698 PRINT "FA -14,50 CF 0 , " ; -LEN( Y 1 abe 1 S ■ .'2* 5/£ ; "D I 0 , 1 . LB " : Yl abe 1 1 . " " 

5500 PRINT "CF C,0 Cl" 

59C4 END IF 
5908 Ipn-0 
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5512 Repeat 1 
5916 Xll-I.E+6 
5920 X u 1 -- 1 . E + E 
5924 Icn»C 
5926 8EEP 

5932 o^-i 

5936 INPUT "WANT TC PLOT OATA FROM A PILE ( 1-Y( OE c AULT ; ,0«N ) ? " ,01 
5940 IP OP « 1 THEN 
5844 BEEP 

5948 INPUT "ENTEF THE NAME OP THE OAT A PILE" , Of l le$( 0 ) 

5952 ASSIGN @File TO Dfile$(0> 

5958 BEEP 
5960 Npairs*20 

5984 INPUT “ENTER THE NUMBER OP X-Y PAIRS STORED ( DEPAULT = 20 ) “ , Npa lr 
5988 BEEP 
5972 Itn-Itn+1 

5978 INPUT "ENTER TU8E NUMBER (1, 2, 3, 4, OR 5>“,Itn 

5980 BEEP 

5984 PRINTER IS 1 

5988 INPUT “UANT DEPAULT SYMBOLS? (YES-0 ( OEP AULT ) , NO* 1 ) “ , Symb 

5992 Syn-Itn+2 

599E IP Synb-e THEN 

6000 60T0 8036 

8004 END IP 



6008 


PRINT USING 


"4X, ““Select a synbol 


6012 


PRINT USING 


"EX , " “ 1 


Star 


2 P 1 u 5 sign""" 


6016 


PRINT USING 


“ 6X , “ “ 3 


Circle 


4 Square""" 


6020 


PRINT USING 


“6X , " "5 


Ronbus “ 


" ■ 


6024 


PRINT USIN6 


"EX, “"6 


Right-s 


ide-up triangle"" 


6028 


PRINT USING 


"EX, "“7 


Up-si de 


-down triangle""" 


6032 


INPUT Sym 








E03E 


PRINTER IS 705 






6040 


PRINT “PU DI 


“ 






8044 


IP Syn-1 THEN PRINT 


“ SM* " 




6048 


IP Syn-2 THEN PRINT 


“ SM+ “ 




6052 


IP Syn-3 THEN PRINT 


"SMo" 




6056 


FOR 1*1 TC Npairs 






6060 


ENTEF ©Pile; 


X y ( * ) 







8064 IF Opo-0 THEN 
8086 Ya * X y ( Itn-1 ) 

6072 Xe*Xy( 1 1 + 1 1 n > 

6077 END IF 

6080 IF Opo-1 THEN 

8084 Ya-Xy(5*!tn ) 

6086 Xe-Xy( 11+Itn) 

6092 ENO IF 

6096 IF Opo-2 THEN 

8100 Ya = X y ( 5*1 tn ) 

6 104 Xa-Xy ( I tn-1 ) 

E 1 08 END IF 

6112 IF Xa <X 1 1 THEN Xll*Xa 

8 1 1 E IF X a > X u 1 THEN Xul=Xe 

8120 IF Iun-1 THEN 

6124 IP Opo .2 THEN Xa=Xa*l.E 

6126 IP Opo>0 THEN >e=>e*.17El 

£132 IF Opc*C THEN Ye*Ye*.3l" 

6136 IP Opc-2 THEN Xa*Xa*.3l7 

6140 END IP 

6144 >-LGT(Xe'Xnin>*£f> 

6146 Y-LGT(Ya/Ywin)*Sfy 
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6152 
6155 
616C 
E1B4 
6 1 BE 
E17: 
E176 
6 1 BC 
61B4 

Biee 

bis: 

6196 

6200 

e:04 

B206 

6212 

6216 

6220 

6224 

6228 

6232 

B236 

B240 

6244 

624B 

6252 

B256 

6260 

6264 

6268 

6272 

6276 

B2B0 

62B4 

6288 

6292 

B296 

6300 



K j * 0 

CALL Symb ( ' , > , S vr. , 1 c 1 , L j ) 

GOTO 6212 

IF Sy r.>2 THEN PRINT *' SH- 
IP Syw'4 THEN PRINT "SR 1. 4,2.4" 

IF lcl-0 THEN 

PRINT -PA* ,X,Y, B " 

E^SE 

PRINT "FA" , X , Y , " P 0 " 

END IF 

IF Syn/3 THEN PRINT “SR 1.2, 1.6" 

IF Syn«4 THEN PRINT ‘1)02 , 4 , 99 , 0 , -8 , -4 , 0 , 0 , £ , 4 , 0 ; " 

IF Sym-5 THEN PRINT "UC3, 0,99, -3, “6, -3, 6, 3, 6, 3, -6. " 
IF Sym*6 THEN PRINT " UC0 , 5 . 3 , 99 , 3 . -8 , -6 , 0 , 3 . 6 , " 

IF 5yn c 7 THEN PRINT " UC0 , “5 . 3 , 99 , “3 , 8 , 6 , C , -3 , -6 ; M 

NEXT 1 

PRINT “PIT 

BEEP 

I lab«i 



INPUT “WANT TO LABEL 7 ( 1 «Y ( DEFAULT ), 0»N )", 1 1 ab 
IF Ilab-I THEN 

PRINT "SP0.SP2" 

BEEP 

IF K 1 ab*0 THEN 
Xlab*E5 
Vlab=85 

INPUT "ENTER INITIAL X,Y LOCATIONS" ,Xlab,Ylab 
Xt t *X 1 ab-5 



Yt t*Y lab + 8 
PRINT "SP 1,1.5- 

PRINT "Sfi.PA" , X 1 1 ,Ytt , " LB Tube X File' 

Yt t»Yt t“3 



PRINT "PA" , X 1 1 , Ytt , "LB No Oil Nane K 

IF Syn«=l THEN PRINT "5n«" 

IF Syn-2 THEN PRINT " 5H+ " 

IF Syn-3 THEN PRINT "SMo" 

Klab=l 
END IF 



E304 k j * 1 

6308 CALL Synb< X lab , V lat , Syn , Ic 1 , k j ) 

6312 PRINT "SF 1 , 1 . 5 ; SM " 

6316 IF Sym<4 THEN PRINT "PR 2,0" 

6220 PRINT “PR 2,0. L8" . I tn ; " " 

6324 BEEP 

6328 INPUT 'ENTER BOP(0»DE cr AULT )" ,Bop 

6332 IF Bop 10 THEN PRINT 'PR 3 , 0 ; LE “ ; Bop . " “ 
6336 IF Bop>9 THEN PRINT “PP 1 . 5 , 0 . LB" ; Bop ; " " 

6340 PRINT "PR :,0.LB";Ofile$<0);"" 

6344 PRINT " S D 0 . S p 1 ; SR 1.5,2" 

E34E Y lab=Yl at-5 

E252 END IF 
6356 EEEP 

6360 ASSIGN ©File TO • 

6384 xii-xii/:.: 

E266 Xul-Xul*:.: 

63“ , 2 I GOTO 8040 
6376 END IF 
6380 PRINT "PU SN" 

£364 BEEP 
6366 Gc_on=l 






) 



) 



) 
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E7S2 INPUT "l )AN T TO PLOT A POlYNOMIA. (1 - r ' DEF AUL T , , P*N : ' . fcc.on 

E 3 9c IP 6c_ on* i THEN 

6400 E£iP 

6404 FRINTEP IS 1 

6406 INPUT "WANT OEFAULT LINE TYPE 7 v YES-0 < OEFAULT >, NG- 1 Ln 

6413 Ipn* I t n 

6416 IF Ln-0 THEN 

6430 SOTO 6448 

6434 ENC IF 

6438 PRINT US INS *4X, '“Select line type 

6433 PRINT USING "6X,“*'0 Solid line”"" 

6436 D R I NT USING "6X,""1 Oashed"“•• 

6440 PRINT USING ■ 6X , " H 2 , , , 5 Longer line - dash 

6444 INPUT Ipn 

6448 PRINTER IS 705 

6453 BEEP 

64SE Ilog-1 

6460 INPUT * SELECT ( 0-LIN , l-LOG( DEFAULT >)", I log 

6464 Iprnt-1 

6468 CALL Poly ( Df l leS( • ) , Npairs , I in ) 

6472 FOR Xx-0 TO Cx STEP Cx/200 

6476 Xa-Xmin* 10 *Xx 

6480 IF Xa<X 1 1 OR Xa>XuI THEN 6572 

6484 Icn-Icn+1 

6488 Pu-0 

6492 IF Ipn-1 THEN Idf-Icn MOD 2 

6496 IF Ipn-2 THEN Idf-Icn MOD 4 

6500 IF Ipn-3 THEN Idf-Icn MOO 8 

6504 IF Ipn-4 THEN Idf-Icn MOO 16 

6508 IF Ipn-5 THEN Idf-Icn MOO 38 

6512 IF Idf-1 THEN Pu-1 

6516 Ya-FNPo 1 y ( Xa ) 

6520 IF Ya<Ynin THEN 6573 

6524 IF Iun-1 THEN 

6528 IF Opo <2 THEN Xa-XaM.8 

6532 IF Opo>0 THEN Ya«Ya*.176l 

6536 IF Opo-0 THEN Ya-Ya*.3l7 

6540 IF 0po=2 THEN Xa-Xa*.3l7 

6544 END IF 

6548 Y-L6T( Ya/Ynin )*Sf y 

6552 X-L6T ( Xa /Xm in )*Sf * 

6556 IF Y<0 THEN Y-0 

6560 IF Y > 1 00 THEN GOTO 6572 

6564 IF Pu-0 THEN PRINT "PA" , X , Y , "PD" 

6568 IF Pu-1 THEN PRINT " PA " , X , Y , " PU " 

6572 NEXT X > 

6576 PRINT "PU M 

6580 ENO IF 
6584 BEEP 

6588 INPUT "WANT TO QUIT ( 1 *Y , B»N( OEFAULT ) ) ” , Iqt 

6592 IF Igt-1 THEN 6600 

6596 GOTO 5916 

6600 PRINT "PU PA 0,0 SP0" 

6604 SUBEND 
6608 1 

6612 ! xxxxxxxxyxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
66 1 6 1 

6620 SUB Synb( X , Y.Svn, Icl ,(■ j ) 

6624 IP Syr - 'Z THEN PFINT "SM" 

6628 IF Sy m < 4 THEN PPINT "SR 1.4, 2. 4" 
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6E32 


Vad-C 




6636 


IF Kj-1 


THEN N ad- . 6 


6640 


IF Icl-0 


THEN 


6644 


PRINT 


H PA* , X , Y4>®d, ” “ 


6646 


ELSE 




6652 


PRINT 


“PA" ,X , Y+Yad, "PD" 


6656 


END IF 




6660 


IF Syr, >3 


THEN PRINT “SR 1 .2 , 1 .6” 


6664 


IF Sym-4 


THEN PRINT "UC2,4,9S,C,-E,-4,0,C,B,4,0 ( " 


6666 


IF Syr-5 


THEN PRINT "UC3,0,9S,-3,-6,-3,E,3,6,3,-6 


6672 


IF Syn-6 


THEN PRINT "UC0, 5. 3, 99, 3, -6, -6, 0,3, 8, " 


6676 


IF Syr-7 


THEN PRINT “ UC0 , -5 . 3 , 99 , -3 , 8 , 6 , 0 , >3 , -8 . H 


6660 


IF Kj-1 


THEN PRINT "SM.PR 0,-.8“ 


6684 


SU8END 




6688' 







6692 'XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
66 96 » 

6700 SUE Fixup 
67041 FILE- FIXUP 
67081 

6712 DIM Enf ( 34 ),Anp{ 11 ) , Uo 1 1 ( 4 ),Ldtc( 4 ) 

6716 BEEP 

6720 INPUT “OLD FILE TD F IXUP M , D2? i le$ 

6724 ASSIGN ©Fi 1 e2 TO D2fileS 
6728 D 1 f i 1 eS* ' TEST " 

6732 CREATE BOAT Dlfile$,60 
6736 ASSIGN ©Filel TD DIfileS 

6740 ENTER © r l 1 e2 : Nrun , Da t e$ , L dt c ( • ) , I t t , 8op , Nht , No t p , Nrt , Corr 
6744 OUTPUT ©F i 1 e 1 ; Nrun , Da t e$ , L dtc ( * ) , It t , Bop , Nht , Na t p , Nr t , Corr 
6748 FOP 1-1 TD Nrun 

6762 ENTEP ©F 1 1 e2 . To 1 d$ , Enf < * ) , Uo 1 1 < • ) , Anp( • ) 

6766 IF 1«1 THEN 6764 

6760 OUTPUT ©Filel.8op,Told$,Emf( * ) , Uo 1 1 ( • ), Amp ( ♦ ) 

6764 NEXT I 

6768 ASSIGN ©File2 TD • 

6772 ASSIGN ©Pile! TD * 

6776 t RENAME “TEST * TO D2_fileS 

67ee 8EEP 2000 , .2 

6784 BEEF 4000, .2 

6768 BEEF 4000, .2 

6792 SU8END 

6796 1 

6800 'XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
6804 ! 

6808 SUB Move 
6£12» FILE NAME MOUE 
68161 

6 620 DIM A(56 > , B ( 66 ),C(66 , D< 66 J , E< 66 ) , F ( 66 ) , G( 66 ) , H( 66 ) , J ( 66 1 , K ( 66 ) , L< 66 ; , M( 66 
) 

682^ DIM N( 66 ) , Enf ( 34 > , Uol t < 2 ) , Amp < 11 ) , Ldt c ( 4 ) 

6826 BEEP 

6632 INPUT “DLD FILE TO MDUE *’ , D2_ f : 1 eJ 
6636 ASSIGN ©F;le2 TO D2_f:Jei 

6640 ENTER ©F 1 1 e2 . Nr jr. , T o 1 ri5 , Ldt c ( ♦ , 1 1 1 .Bop , Nhi , Netp , , Ccrr 

6844 FDP 1-1 TO Nrur 

6548 ENTEF ©File!. Told! 

6852 ENTEF ©F 1 1 e2 . A( 1 ) , B< I ) , CM ' , D< I > , E Cl j , F« I \ G( I • . Hi I > , J ' I , J - 1 > , L • I ■ , M< 

1 ) ,N< I ) 

6e5E ENTEF ©F 1 1 el ; Enf ( * : , Uo 1 1 ' ♦ ; , Anp ( « ) 

6650 NEXT 1 
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ES64 TC * 

6656 EEE C 

66": I NPUT " S H I r T DISK AND HIT CONT I NUE “ , DK 

6676 EEE C 

6 88C INPUT “INPUT BDAT SIZE' , S l z e 
6864 CREATE 6DAT D2_f i let , S: re 
6 88E ASSIGN e-ilel TO D2_file$ 

68SC OUTPUT @F i le 1 ; Nr U n , Da t eS , Ldt C < • ) . 1 1 1 , Bop , Nht , No t p , Nr t , Cor r 
6896 FOR I«1 TC Nr u r, 

6900 OUTPUT i I e 1 ;TcIdS 

6 904 OUTPUT @Filel.A(I),B(I) t C<l),D(I),E<I),F(I>,G(I>,H(I),J(l7,K(I),L<l>,M 

< I ) ,N( I ) 

6908 DUTPUT GFi lei : Emf ( « > , Uo 1 1 < • > , Amp ( • } 

6912 NEXT I 

6916 ASSIGN @F 1 1 e 1 TD * 

6920' RENAME "TEST" TO D2_fiIeS 

6924 BEEP 2000, .2 

6926 BEEF 4000, .2 

6952 BEEP 4000, .2 

6936 PRINT "DATA FILE MDUED" 

6940 SU8END 
69441 

6948 I xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
6952 r 

6956 SUB Purg 
6960 BEEP 

6964 INPUT "ENTER FILE NAME TD BE DELETED ", F l leS 
6968 PURGE Files 
6972 GDTD 6960 
6976 SUBEND 
6 980 1 

6984 l xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
6988' 

6992 SUB Comb 

69961 FILE NAME COMB 

70001 

7004 DIM Emf<34 ),Uolt<2 \Amp( 11 >,Ldtc< 4 ) 

7008 BEEF 

7012 INPUT "OLD FILE TG FIXUP" , D2_f i leS 
7016 ASSIGN @F i le2 TO D2_fiieS 
7020 D l_f l le$* " TEST " 

7024 CREATE BDAT Dl_fileS,30 
7028 ASSI6N SPilel TO Dl_fileS 

7032 ENTER @F ; 1 e2 . Nrun , DateS, Ldtc(* ), 1 1 1 ,Bop,Nht , Na t p , Nr t , Corr 

7036 IF K*0 THEN DUTPUT @F i Ie 1 ; Nrun , Da t e$ , Ldtc < * ) , 1 1 1 , Bop , Nht , Na tp , Nr t , Corr 

7040 FOR 1*1 TO N r un 

7044 ENTER @F i le2 ; Bop , To 1 dS , Er f ( • 5 , Uo 1 1 ( ♦ ) , Arp ( * ) 

7046 OUTPUT £F i le 1 ; Bop , To! dS , Emf ( • ) , Uol t < * > , Amp ( • ) 

7052 NEXT I 

705E ASSIGN @-:Ie2 TO * 

7060' RENAME “TEST" TO D2_f:leS 
7064 BEEP 4000, .2 
7068 BEEP 
70.72 OK e = 1 

707E INPUT “UAN T TO ADD ANOTHER FILE (1*Y,0*N< default ) ,Ot s 
7080 IF 01 ?=1 THEN 
7064 K-I 
70B8 BEEP 

7092 INPUT "G1UE NEU FILE NAME " , Nf i 1 eS 
7096 ASSIGN §Pile2 TO NfileS 
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7100 GOTO 7030 
7104 END IF 

710£ ASSIGN PFileZ TO * 

71 IS SUBEND 
7116' 

7i:e ' xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

7134 i 

7128 SUE Readplot 

7132 DIM Qdp ( 5 > , Ht ube ( 5 ) , The t ab( 5 ) 

7136 PRINTER IS 701 

7140 INPUT ‘ENTER FILE NAME * , F i 1 e* 

7144 INPUT “ENTER THE NUMBER OF DATA PAIRS”, Nr U n 
7146 ASSIGN CFilel TO File$ 

7152 FOR 1-1 TO Nr U n 

7156 ENTER 5F l le 1 : Odp ( • ) , Ht ube < • ) , The t ab< * ) 

7160 PRINT Odp ( • ) 

7164 PRINT 

7166 PRINT Htube ( * ) 

7172 PRINT 

7176 PRINT Thet ab< ♦ ) 

7180 PRINT 

7184 PRINT 

7188 NEXT 1 
7192 SUBEND 
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